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Abstract 

TLe  parpos^  of  this  research  was  to  devdop  a  procedure  tot  estimating  the 
average  histatio  radar  dntter  power  from  a  topographic  map  and  using  rough  surface 
scattering  theory,  first,  the  terrain  in  the  radar  target  area  was  divided  into  sub 
regions  having  nearly  homograeous  diaracteristics.  A  joint  Gaussian  distribution  for 
the  surface  hdghta  was  fitted  from  the  information  contained  in  the  topographic  map 
contour  fines  fr>r  each  homogeneous  area.  Bcom  these  distributions,  the  normalized 
radar  cross  section  fcr  each  terrain  area  was  determined.  The  clutter  power  for  each 
range  cell  was  then  determined  via  a  modified  radar  range  equation. 

The  procedure  was  implemented  for  three  target  areas  of  an  existing  ground 
based  Hstatic  system  near  Hanscom  AFB,  Massachusetts  using  a  U.  S.  Geological 
Service  topographic  map  of  the  area.  Scattered  power  from  very  rough  surfaces  was 
estimated.  However,  the  accuracy  of  the  estimates  could  not  be  verified  due  to  a  lack 
of  actual  data  for  comparison.  It  was  also  (fetermined  that  topographic  map  data 
may  not  be  sofiicient  to  determine  scattered  power  from  slightly  rough  surfaces. 
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DETERMINATION  OF  BIST)\TIC  RADAR  CLUTTER  POWER 
FROM  TOPOGRAPHIC  DATA 


L  Introduction 

1.1  Background. 

The  operating  principles  of  radar  are  relativelj  simple.  The  process  initiates 
when  a  transmitter  emits  an  identifiable  dectromagnetic  wave  into  the  atmosphere. 
'When  the  wave  strikes  an  object,  electromagnetic  energy  scatters  in  all  directions.  A 
receiver  dstects  the  presence  of  a  target  when  the  reflected  energy  strikes  the  receive 
anteima.  Given  knowledge  of  the  exact  time  of  transmission  of  the  original  wave 
and  of  the  vdodty  of  wave  propagation,  the  receiver  can  calculate  the  approximate 
range  between  the  target  and  the  radar  system,  the  energy  from  the  transmitter  is 
concentrated  in  a  narrow  beam,  the  bearing  and  elevation  of  the  target  can  also  be 
determined.  With  these  capabilities,  it  is  easy  to  see  how  radar  has  evolved  into  a 
necessary  piece  of  military  hardware.  It  finds  application  on  every  military  aircraft 
in  the  inventory  today  and  its  capabilities  often  play  a  vital  role  in  the  success  or 
frulure  of  a  mission. 

Problems  arise  when  the  reflected  energy  is  not  due  soldy  from  a  desired  target. 
'When  the  transmitting  anteima  directs  energy  at  low  grazing  angles,  objects  on 
the  ground  reflect  energy  to  the  receiver.  This  phenomenon  is  known  as  dutter. 
The  amount  of  dutter  energy  returned  to  the  radar  receiver  is  proportional  to  the 
reflectivity  of  the  terrain  and  other  &ctors.  'When  more  energy  is  reflected  from 
terrain  flsatures  than  from  a  desired  target  the  capabilities  of  a  radar  are  reduced. 


All  military  forces  attempt  to  tise  dutter  to  their  advantage.  Low  flying  air¬ 
craft  and  missiles  are  difficult  to  detect  in  the  presence  of  dutter.  The  dutter  energy 
obscures  target  reflectioxis  and  provides  an  attacking  force  with  the  advantage  of  less 
time  of  exposure  to  hostile  countermeasures.  The  detection  of  low  flying  objects  is 
compounded  by  the  use  of  stealth  technology,  which  reduces  the  amount  of  energy 
reflected  in  the  backscatter  direction.  Monostatic  radar  ^sterns  use  collocated  trans¬ 
mitters  and  receivers.  The  transmitter  and  receiver  often  share  a  common  antenna 
and  the  system  only  detects  backscattered  target  energy.  With  less  dectromagnetic 
energy  reflected  toward  a  receiver  by  using  stealth  technology,  the  task  of  hiding  an 
incoming  vehide  in  ground  dutter  becones  easier. 

This  work  is  sponsored  by  the  Applied  Electromagnetics  Division  f  the  Rome 
Air  Devdopment  Center  (RADC).  RADC  has  proposed  using  bistatic  radar  ^sterns 
with  spedal  oonflgurations  as  a  techniqiie  to  enhance  the  detection  of  low  observable 
(stealthy)  objects.  Bistatic  radar  systems  use  transmitters  and  receivers  which  are 

I 

not  collocated  (18:1)  and  do  not  rdy  oh  backscattered  energy  for  proper  operation. 

i 

RADC  proposes  that,  for  a  given  terrain  and  bistatic  angle,  there  will  be  a 

transmitter  to  recdver  antenna  polarization  relation  which  will  result  in  a  minimum 

! 

amount  of  dutter  return  from  the  terrain  in  question.  Use  of  the  process  to  enhance 
the  detection  of  bw  observable  objects  requires  apriori  knowledge  of  the  required 
antenna  polarization  relationship  which|  is  obtained  from  the  reflective  characteristics 
of  the  terrain.  To  obtain  optimum  perfonnance  of  the  Ihstatic  oonflguration  proposed 
RADC,  a  thorough  knowledge  of  the  dutter-to-noise  power  profiles  of  bistatic 
radars  operating  at  bw  grazing  ang'  '  is  needed. 

For  the  foregoing  reasons  there  is  considerable  interest  in  predicting  the  terrain 
reflectivity  seen  by  Ihstatic  radar  systems  at  bw  grazing  angles.  The  dutter-to-noise 
power  levds  at  a  radar  receiver  antenna  are  directly  related  to  the  reflectivity  of  the 
terrain  in  the  target  area.  Several  analytical  and  empirical  models  (13:671-773),(9:74- 
134), (3:70- 119)  have  been  proposed  to  calculate  the  terrain  refl< activity  from  known 
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t c'poloi'i.-nj  fratuic*. 

Tlio  amount  of  clutter  power  a  bistatic  radar  detects  depends  on  the  radar 
i  arameters,  the  gfx)metric  bcations  of  the  transmitter,  receiver,  and  target  area  as 
V  -  il  .xs  t'-.e  I'liysical  characteristics  of  the  terrain. 

Tiie  radar  parameters  include  the  frequency  of  operation  of  the  transmitter, 
the  ttansmit  and  receive  antenna  polarizations,  the  transmit  and  receive  antenna 
Ix  ainwidths,  and  the  transmitter  pulse  width.  All  of  these  parameters  are  control- 
Lil  le  by  the  radar  design  engineer. 

The  geometric  considerations  involve  the  physical  locations  of  the  transmitter 
and  the  receiver  in  relation  to  the  target  area.  The  most  prominent  of  these  factors 

are; 

•  The  distance  between  the  transmitter  and  receiver. 

•  The  distance  from  the  transmitter  to  the  target  area. 

•  The  distance  from  the  receiver  to  the  target  area. 

•  Tlic  bistatic  angle  formed  by  the  transmitter,  receiver,  and  target. 

•  The  height  of  the  transmitting  and  receiving  antennas. 

•  The  grazing  angle  of  the  transmitting  and  receiving  antennas. 

Tlic  physical  characteristics  of  the  terrain  can  not  be  controlled  and  must  be 
measured  or  estimated  hr  each  particular  ^)plication.  The  terrain  factors  which 
alTect  the  amount  of  reflected  duttcr  power  are: 

•  The  physical  ciiaracteristics  of  the  terrain. 

1.  Flat  land,  rolling  hills,  mountains,  bodies  of  water. 

2.  Tlie  amount  of  moisture  in  the  soil. 

3.  Tlie  dielectric  constant  of  the  terrain. 
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•  The  population  of  the  terrain. 


1.  The  spediic  type  of  vegetation,  trees,  or  buildings  covering  the  terrain. 

2.  The  relative  heights,  thicknesses,  and  densities  of  the  population. 

3.  The  dielectric  constants  of  the  population  items. 

•  The  weather  conditions  in  the  target  area. 

All  of  these  factors  affect  the  amount  of  clutter  power  received  by  a  bistatic  radar 
system. 

1.2  Problem  Statement. 

No  accurate  procedure  exists  for  determining  the  dutter-to>noise  power  profile 
of  a  ground  based  bistatic  radar  system  operating  at  low  grazing  angles  and  large 
Bistatic  angles  from  knowledge  of  the  terrain  features  of  the  radar  target  area. 

l.S  Research  Objective. 

The  purpose  of  this  research  effort  is  to  develop  an  accurate  procedure  fer 
predicting  the  clutter-to-noise  power  profile  for  a  ground  based  bistatic  radar  system 
from  knowledge  of  the  target  area  topography. 

1.4  The  Bistatic  System. 

The  transmitter  fer  the  bistatic  radar  system  used  fer  this  thesis  was  an  S  band 
weather  radar  located  near  Stow,  Massachusetts.  The  receiver  was  an  RADC  asset 
located  atop  Prospect  lEll  in  Waltham,  Massachusetts.  The  proposed  target  area 
was  the  terrain  southwest  of  Hanscom  Air  Fbrce  Base.  Stow  is  located  approximately 
20  kilometers  (Km)  west  of  Waltham.  RADC  originally  requested  dutter-to*noise 
power  ratios  fer  five  range  bins  at  three  different  bistatic  angles  between  120*  and 
160*.  The  clutter  powers  were  to  include  the  effects  of  the  power  associated  with  the 
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main  lobes  and  significant  side  lobes  of  both  antennas.  The  researcher  was  at  liberty 
to  select  the  appropriate  bistatic  angles  and  range  bins.  As  the  project  progressed, 
three  specific  sites  were  selected  by  RADC. 

The  transmitter  was  operated  at  a  frequency  of  3  gigahertz  (GHz),  had  a 
peak  power  of  800  kilowatts  (KW),  a  pulse  width  of  1  microsecond  (ps),  and  a 
pulse  repetition  frequency  (PIIF)  of  1  kilohertz  (KHz).  The  transmitting  antenna 
was  a  28  foot  parabolic  reflector  and  could  be  operated  with  vertical  or  horizontal 
polarization. 

The  range  gated  receiver  was  equipped  with  a  20  dB  standard  gain  horn  an¬ 
tenna  and  could  operate  with  cither  vertical  or  horizontal  polarization.  The  term 
standard  gain  horn  was  taken  to  mean  a  p3rramidal  horn  antenna  designed  for  opti¬ 
mum  gain  in  both  the  E  plane  and  H  plane. 

1.5  Assumptions. 

The  radar  range  cells  were  assumed  to  be  in  the  far  field  of  both  the  transmit¬ 
ting  antenna  and  receiving  antenna.  This  implied  the  vravefronts  incident  upon  the 
terrain  and  the  receiving  antenna  aperture  were  planar.  This  vras  a  valid  assump¬ 
tion  given  the  geometry  of  the  test  area.  Skolnik  (15:229)  describes  the  far  field  as 
a  distance  greater  than  Rp  where  Rp  is  defined  as: 

Rp  =  D^X  (1.1) 


Where: 

D  =  the  largest  dimension  of  the  antenna  aperture  in  meters  (m) 

A  =  the  wavelength  o^  the  radar  in  meters 
The  wavelength,  A,  was  determined  from  (11:1) 

X-c/f  (1.2) 


1-5 


Where: 


c  =  the  velocity  of  wave  propagation,  c  =  3X10*  meters/second  (m/s) 

/  =  the  operating  frequency  of  the  radar  in  ^rtz  (Hz) 

Using  an  (derating  frequency  of  3.0  GHz  in  Eq  1.2  yielded  A  =  10  centimeters 
(cm).  The  largest  dimension  of  the  receiver  antenna.  ^>ertnre  was  45  cm  and  the 
transmitter  antenna  had  an  diameter  of  28  feet  or  8.5344  m.  Eq  1.1  was  used  to 
calculate  Rp  as  2.025  m  fer  the  receiver  and  728.36  m  fer  the  transmitter.  Since  the 
range  from  either  the  transmitter  or  the  receiver  (20  Km)  was  much  greater  than 
Rp,  the  target  area  was  in  the  fer  field  and  the  planar  wavefront  assumption  was 
valid. 

The  procedure  developed  during  this  research  was  limited  to  stationary  ground 
based  bistatic  systems  operating  at  low  grazing  angles  and  all  bistatic  angles  except 
the  ferward  scatter  and  back  scatter  directions.  The  radar  receiver  was  assumed  to 
be  operating  without  using  pulse  integration. 

1,6  Scope. 

This  thesis  was  limited  to  the  development  of  a  process  fer  the  determination 
of  the  dutter-to-noise  power  profile  fer  a  generic  ground  based  bistatic  system.  The 
process  was  then  applied  to  five  range  cells  fer  three  bistatic  angles  of  the  specific 
Ustatic  radar  ^stem  used  by  RADC.  The  range  cells  selected  fer  analysis  represent 
only  a  small  sample  of  the  total  target  area.  The  dutter-to-noise  power  profile 

concept  could  easily  be  extended  to  an  entire  target  area. 

1 

The  effects  on  the  dutter-to-noise  powp  profiles  of  extreme  weather  conditions, 
variations  in  soil  moisture  content, and  variations  in  target  area  population  were  not 
considered  in  this  research  effort. 
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1.7  Approach. 

The  lesearch  was  accomplished  in  several  steps.  Barxick  (13:671-729)  details 
processes  ibr  calcolating  scattered  power  from  two  dasses  of  rough  surfaces.  His 
description  encompasses  coherent  and  incoherent  scattering  from  slightly  rough  sur- 
&ces  and  incoherent  scattering  from  very  rough  surfaces.  Both  cases  treat  the  height 
of  the  surface  as  a  random  variable  and  require  knowledge  of  the  distribution  of  the 
surface  heights.  Therefore  each  target  site  was  subdivided  and  analyzed  to  obtain 
the  necessary  statistical  information. 

Topographic  data  fi3r  the  sdected  target  sites  was  digitized  and  entered  in  a 
data  base.  The  data  was  analyzed  to  develop  statistical  representations  of  the  terrain 
characteristics.  Ebdsting  rough  suxftu:e  scattering  procedures  were  used  to  predict 
the  dutter-to-noise  power  profile  of  three  range  cells  at  each  site.  The  process  was 
initiated  by  determining  the  footprint  of  the  bistatic  radar. 

1.7.1  Radar  Antenna  Footprints.  The  radar  footprint  is  the  terrain  area  de¬ 
lineated  by  the  intersection  of  the  projected  transmit  and  receive  antenna  patterns 
and  the  concentric  dhpses  which  define  the  isorange  contours  of  a  bistatic  radar 
qrstem.  The  size  and  shape  cS  the  individual  range  cells  were  determined  by  several 
factors: 

1.  The  radiation  patterns  of  both  antennas. 

2.  The  relative  heights  of  the  antennas  and  the  target  area. 

3.  The  bistatic  angle  formed  by  the  transmitter,  receiver,  and  target  area. 

4.  The  distance  from  transmitter  to  target  and  receiver  to  target. 

5.  The  range  gating  interval  used  in  the  receiver. 

1.7.2  Terrain  Maps.  Topographic  maps  &r  the  target  area  were  obtained 
from  the  United  States  Geological  Sendee  (USGS)  to  correlate  the  physical  diarac- 
texistics  of  the  terrain  to  the  selected  range  cdls.  The  maps  were  scaled  to  1/250000 


\ 
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and  allowed  analysis  of  the  grotind  cover  in  the  area  of  interest.  The  three  range  cells 
selected  ibr  evaluation  at  each  l^tatic  angle  were  subdivided  into  smaller  blocks  £)r 
statistical  analysis  to  determine  the  mean,  variance,  and  correlation  length  of  the 
terrain  height  distributions  covered  by  the  selected  range  cells. 

The  information  from  the  topographic  maps  fcr  the  specific  blocks  of  each 
range  cell  was  transferred  to  a  memory  array  where  the  data  was  used  to  determine 
the  large  scale  of  roughness  parameters.  Specifically,  the  data  represents  the  relative 
position  (x  -  y  coordinates),  height,  didectric  constant,  and  type  of  the  terrain 
inhabiting  the  range  cell  blocks. 

Though  the  USGS  maps  are  detailed,  they  lack  sufficient  information  to  de¬ 
scribe  the  range  cell  blocks  completely.  The  maps  provide  detailed  information  on 
the  locations  and  relative  sizes  of  buildings,  roads,  streams,  and  major  bodies  of 
water.  The  elevation  data  is  provided  in  the  fbim  of  three  meter  contour  fines  and 
is  accurate  enough  to  determine  the  necessary  large  scale  of  roughness  parameters 
required  for  the  RADC  application. 

The  land  cover  on  the  USGS  maps  was  identified  as  frdling  into  one  of  ax- 
teen  broad  categories.  Examples  of  these  categories  are  woodland,  scrub,  marsh, 
swamp,  mangrove,  orchard,  and  vineyard.  Though  informative,  the  information  was 
not  detailed  enough  to  determine  the  small  scale  of  roughness  parameters  for  the 
specific  RADC  ^plication.  If  the  entire  range  cdl  in  question  had  a  homogeneous 
population,  eg.  all  of  the  blocks  in  range  cell  one  were  entirdy  covered  by  pine 
trees,  additional  information  would  still  be  required.  Hie  additional  information 
would  have  to  be  detailed  enough  to  identify  individual  tree  heights  so  thsd  a  height 
distribution  could  be  developed.  T^^e  mean  and  variance  of  the  height  distribution 
would  then  be  used  to  calculate  the  incoherent  power  scattered  by  the  small  scale  of 
roughness.  Information  is  also  needed  to  determine  the  relative  permittivity  of  the 
target  area  surface  and  any  vegetation  which  may  inhabit  surface. 

Unfortunately  nature  is  not  this  cooperative.  It  is  unlikely  that  the  entire 
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nmge  cell  will  have  a  homogeneous  small  scale  of  roughness.  The  extreme  case 
would  inTolve  a  different  small  scale  of  roughness  for  each  range  cell  block.  In  this 
instance  a  separate  height  distribution  would  have  to  be  developed  for  each  range  cell 
block.  This  process  could  rapidlj  become  computationally  cumbersome  if  developing 
a  dutter-to-noise  power  map  for  the  entire  coverage  area  of  a  spediic  bistatic  radar 
^stem. 

1.7.S  Clutter-to~Noise  Power  Profile  Maps.  To  simplify  the  calculations  re¬ 
quired,  a'eas  of  homogeneous  surface  diaracteristics  were  selected  at  each  target 
area.  The  height  distributions  in  these  areas  were  assumed  to  be  Gaussian  in  na¬ 
ture.  The  mean,  variance,  and  correlation  distance  were  extracted  from  the  data 
available  on  the  topographic  maps.  The  characteristics  of  the  identified  Gaussian 
PDF  were  used  to  implement  Barrick’s  rough  surface  scattering  frrmulas. 

Applying  Barrick’s  formulas  to  each  range  cell  diaracterization  area  provided 
a  value  for  the  normalized  radar  cross  section  (NRCS)  of  the  block.  The  area  of 
the  selected  range  cells  was  used  to  determine  an  average  value  of  RCS  which  in 
turn  vras  used  in  a  variation  of  the  Ustatic  radar  range  equation  to  determine  the 
scattered  power  at  the  receiver,  ha  the  cases  where  the  terrain  exhibits  both  scales 
of  roughness  over  a  single  range  cell,  the  scattered  power  for  both  cases  could  be 
calculated  and  added  vectorially. 

The  implementation  of  the  scattering  formulas  was  essentially  the  same.  The 
main  difference  in  the  two  cases  (slightly  rough  surfaces  and  very  rough  surfaces) 
was  in  how  they  were  developed  and  their  surface  height  variations  relative  to  the 
wavelength  of  the  excitation  source.  The  very  rough  surface  formulas  were  developed 
using  physical  optics  principles.  The  slightly  rough  surface  formulas  were  developed 
using  a  pertiubation  technique. 

Hotting  the  values  of  ducter-to-noise  power  for  each  cdl  relative  to  the  cell’s 
physical  position  yields  a  dutter-to-noise  power  profile  map  for  the  bbtatic  radar. 
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The  concept  is  easily  extended  to  include  the  entire  radar  footprint. 

1.7.4  Performance  Evaluation,  The  results  of  this  research  were  to  be  verified 
with  the  equipment  already  in  place  at  the  RADC  Instatic  radar  sites.  The  actual 
data  was  unavailable  at  the  time  of  publication. 

1.8  Overview  of  the  Thesis. 

Chapter  11  presents  the  theory  behind  the  procedures  pertaining  to  the  deter¬ 
mination  of  clutter  reflectivity.  Chapter  III  presents  the  results  of  the  research  effort 
and  Chapter  IV  presents  conclusions  firom  the  thesis  and  malces  recommendations 
for  further  research. 
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IL  Theory 


2.1  Overview. 

This  chapter  describes  the  procedures  used  to  calculate  the  dutter-to-noise 
power  ratios  £}r  each  range  Hn  and  provides  the  theory  behind  the  concepts  employed 
here. 

2.2  Clutter-To-Noise  Power  Ratio. 

la  most  radar  applications,  the  dutter  power  received  by  a  qrstem  is  unwanted 
and  often  treated  as  noise.  In  this  application,  the  dutter  power  (Pe)  was  treated  as 
the  return  from  a  desired  target.  The  dutter  power  frr  a  bistatic  tystem  given  by 
Skolnik  (15:557). 


PtGtGrX^ffi 


Where: 

Pr  =dntter  power  received 

Gt  —transmitting  anterma  gain  in  the  direction  of  the  dutter 
Gr  =reodving  antenna  gain  in  the  direction  of  the  dutter 
CTi  i=bistatic  radar  cross  section  of  the  dutter 
Dt  j=transmitter  to  target  distance 
Dr  ==receiver  to  target  distance 

Lp{ti  =propagation  loss  over  the  transmitter  to  target  path 
Xp(r)l=propagation  loss  over  the  receiver  to  target  path 
Lg  =system  losses 


(2.1) 
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A  majority  of  the  parameters  in  Equation  2.1  are  known  or  can  be  caicuiated 
given  enough  information  about  the  transmitting  hardware,  receiving  hardware,  and 
physical  layout  of  the  bistatic  system.  The  only  factor  intimately  associated  with 
the  terrmn  characteristics  is  trj  which  represents  the  bistatic  RCS  of  the  target. 
For  many  simple  geometric  shapes  the  RCS  is  known  or  can  be  calculated  with 
reasonable  accuracy.  The  problem  at  hand  is  that  the  target  of  interest  is  terrrdn  of 
varying  degrees  of  roughness  populated  by  objects  of  random  shapes,  sizes,  heights 
and  densities. 

The  randomness  of  the  problem  requires  a  probabilistic  approach  to  determin¬ 
ing  a  value  for  aj  in  realistic  terrains  encountered  by  a  bistatic  radar  system.  The 
value  obtained  for  o-j  is  directly  related  to  the  physical  dimensions  of  the  clutter 
patch  illuminated  by  the  transmitted  signal  and  simultaneously  visible  to  the  re¬ 
ceiving  antenna’s  aperture.  Once  the  visible  footprint  is  identified,  other  factors  can 
influence  the  amount  of  clutter  power  received  from  the  footprint. 

2.2,1  Factors  Affecting  the  Clutter  Power.  The  other  factors  which  affect  the 
amount  of  clutter  power  received  include  the  type  of  land  cover  found  in  the  footprint 
area,  the  density  of  the  land  rover,  the  variation  in  the  surface  heights  and  slopes  of 
the  land  as  well  as  the  variation  in  the  heights  and  rl  'pes  of  the  objects  populating 
the  surface. 

In  addition  to  the  physical  locations  of  objects  in  the  target  are.a,  one  must 
also  consider  the  dectrical  properties  of  these  objects.  The  dielectric  constants  of 
both  the  terrain  surface  in  a  dutter  cell  and  of  the  objects  populating  the  cell  affect 
the  average  amount  of  scattered  |x>wer  reaching  the  receiver  aperture. 

The  average  dutter  power  received  will  also  depend  upon  the  incident  angle  of 
the  illuminating  wave  (dj),  the  angle  of  departure  of  the  scattered  wave  {0,),  and  the 
bistatic  angle  (/?)  formed  by  the  transmitter,  receiver,  and  target  area.  The  angles 
6i,  0,,  and  <f),  are  defined  at  the  point  of  incidence  as  shown  in  Figure  2.1  and  the 
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Figure  2.1.  Bistatic  scattering  geometry, 
bistatic  angle,  /9,  is  defined  as  t  ~ 

Still  other  &ctors  can  affect  the  amonnt  of  scattered  energy  received  from  the 
datter  cell.  The  homogeneity  of  the  surface  composition  in  the  dutter  cell  can  affect 
the  agnal  contribution  from  the  slightly  rough  and  very  rough  surface  contributions. 
The  size  of  the  dutter  cell  blocks  can  have  the  same  effect  on  the  signal  levels  as 
the  homogeneity  of  the  surface.  Surfaces  which  have  small  homogeneous  sections 
of  many  different  types  of  surface  present  a  problem  when  trying  to  determine  the 
dominant  material  in  a  dutter  cell  block.  The  same  situation  is  possible  if  the  size 
of  the  dutter  cell  blocks  are  too  large. 

Weather  conditions  can  diange  the  reflectivity  of  an  area  in  a  matter  of  minutes. 
'Vl^dy  conditions  can  cause  the  RMS  height  of  crops  to  vary  drastically  or  cause 
portions  <ff  trees  to  move  in  and  out  of  a  particular  range  cell  block.  Rain  can  affect 
most  surface  areas  by  dianging  the  surface  moisture  content.  Snow  and  ice  cover 
can  totally  diange  the  reflection  characteristics  of  a  surfrure  area. 

Some  factors  cause  seasonal  change.  Areas  covered  by  vegetation  will  have  an 
RMS  height  which  will  change  as  the  vegetation  grows  over  the  summer  or  is  har> 
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vested  in  the  iall.  Similaxly,  forested  areas  of  deddaous  trees  will  have  a  reflectivity 
which  will  change  seasonally  as  leaves  appear  in  ^ring,  grow  throu^ont  the  summer 
and  drop  off  in  autumn. 

The  previous  discussion  was  intended  to  point  out  that  there  is  no  single  solu¬ 
tion  to  the  problem  of  predicting  the  bistatic  reflectivity  of  a  given  area.  So  many 
factors  can  affect  the  results  that  it  is  only  fisasible  to  estimate  the  value  of  clutter 
power  for  a  given  set  of  conditions.  The  estimate  must  be  revised  as  the  conditions 
diange.  The  first  obstacle  to  overcome  in  providing  this  estimate  is  identifying  the 
clutter  cell  location. 

2.3  Determining  the  Antenna  Footprint. 

The  terrain  area  which  scatters  electromagnetic  energy  that  teaches  the  receiver 
during  any  particular  range  gate  period  is  of  primary  interest.  Determination  of 
these  terrain  sections  is  the  foundation  for  all  subsequent  scattering  calculations. 
Each  dutter  cell  is  delimited  by  two  fsatures,  the  bistatic  eqnirange  contour  lines 
for  a  specific  receiver  range  gate  period  and  the  intersection  of  the  transmitter  and 
receiver  antenna  patterns.  The  equirange  contour  lines  establish  the  areas  of  possible 
energy  return  during  a  range  gate  period.  The  transmitter  and  receiver  antenna 
patterns  are  then  used  to  identify  subsections  of  the  range  odl  vdiich  nudce  significant 
contributions  to  the  received  dutter  power. 

2.3.1  The  Equirange  Contour  Lines.  Figure  2.2  represents  a  two  dimensional 
view  of  the  bistatic  geometry  from  which  it  is  desired  to  determine  the  dimensions 
of  the  individual  range  cells.  In  this  diagram  the  transmitter,  receiver,  and  target 
are  in  the  same  plane  and  variations  in  surface  height  are  being  neglected  for  the 
moment.  This  is  the  same  situation  one  is  &ced  with  when  extracting  data  from  a 
topographic  map. 

The  angle  P  is  the  bistatic  angle.  The  segment  Db  is  the  basdine  distance  in  the 
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figure  2.2.  Twc  dimensional  bistatic  geometry. 


xy-plaae  from  transmitter  to  leceiTer  and  remains  constant  for  a  given  ground  t)ased 
Instatic  system.  The  segment  Dr  is  the  distance  in  the  sgr-plane  from  transmitter 
to  a  selected  range  cell.  The  segment  Dr  is  the  distance  in  the  :qr-plane  from  the 
selected  range  cell  to  the  receiver.  Both  segments,  Dr  and  Dj,  will  vary  in  length 
with  different  range  cells. 

The  amount  of  dutter  power  detected  by  the  receiver  from  eadi  range  cell  will 
be  equal  to  the  average  of  the  energy  intercepted  ty  the  teedver  during  a  single 
range  gate  period.  Therefore  the  width  of  the  range  odl  is  related  to  the  range  gate 
interval  of  the  receiver.  The  number  cf  range  cells  available  is  controlled  by  the 
transmitter  PRF.  The  easiest  way  to  visualize  this  is  to  examine  a  dining  diagram 
for  the  system. 

Figure  2.3.1  is  a  timing  diagram  for  the  general  ground  based  bistatic  system. 
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Figure  2.3.  Ground  based  bistatic  system  timing  diagram. 

The  timing  diagram  is  referenced  to  the  receiver.  At  time  To,  energy  travding  along 
a  direct  path,  the  distance  Rb,  reaches  the  receiver.  Where  Rb  is  the  Ene  of  sight 
(LOS)  or  slant  range  distance  from  the  transmitter  to  the  receiver.  Li  Hgure  2.3. 1| 
tr  represents  the  range  gate  interval  of  the  receiver.  The  energy  £>r  range  cell  one 
strikes  the  receiver  between  time  To  and  To  +  tr.  Similarly  the  energy  for  range  cell 
two  strikes  the  tecdver  antenna  between  To  +  tr  and  To  +  2tr.  The  total  number 
of  range  cells,  n,  is  determined  by  the  PBP  of  the  transmitter.  Li  Equation  2.2| 
Tf  =  has  been  used. 


”  (FRF){tr)-7r 

I 

I 

The  shape  of  a  dutter  cell  for  a  lustatic  system  is  more  complicated  thah 
its  monostatic  counterpart.  The  most  direct  method  hr  determining  the  shape  of 
a  dutter  cell  is  to  examine  the  equirange  contour  lines  for  a  bistatic  system.  For 
any  arbitrary  equirange  contour  line,  energy  reaches  the  receiver  after  a  time  of 
To  +  titr  seconds.  The  distance  traveled  by  this  energy  is  equal  to  the  sum  of  the 
transmitter  to  target  and  target  to  receiver  distances.  An  equivalent  LOS  distance 
can  be  determined  from  the  wave  speed  and  time  of  travd.  At  the  time  of  the 
range  gate,  the  wave  has  been  traveling  for  a  period  of  To  +  titr  seconds..  Equating 
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’  i.Muc  2.1.  l.ition.shij)  Ix^twwii  slant  ran>;cs  and  xy-plane  distances. 
l!;(  Cu'.ts  yields: 

I'f  +  Rn  =  c(T’o  +  nr/j)  =  /?»  +  ncr/i  (2-3) 

'i  'e'  li  'lit  hand  side  of  Iv<ination  2.3  Is  a  constant  value  for  any  ;;ivcn  value  of 
ri  .cid  T/j,  tlie  leeciver  ran-'e  j',at<;  interval.  If  the  tran.smittcr  is  assumed  to  be  an 
i  eltojdc  r.idi.ilor,  r-ner  'y  i.s  radiatfxl  in  all  directions.  Dy  pl()ttin(?  the  target  points 
f  r  (dl  dirrrtions  radi.dly  fro:;i  the  transmitter  and  iletermining  the  values  for  /?r 
•■'.iid  /I’/j  whii  ’i  rati  fy  l/pt  ition  2.3,  tlie.s<r  r<]uirange  contour  lines  will  be  elliptical 
in  eh  ape.  lio’veyer  /I’j-  a;id  ;i’;j  are  the  .slant  ranges  from  the  transmitter  to  target 
.■:;d  lecf  ivfT  to  t.n  p  t  resjH-ctively  and  can  not  1)0  directly  [dotted  on  a  topographic 

!•;  p. 

I  i  ’nir'  11.4  !!!u'.tr.ati  :i  the  ii'I.itionsliii)  iM'twix'n  the  LOS  .sl.ant  ranges  Rt, 
!'■))  .'iiid  Ih*  ir  coricNpoiniiiig  xy-j'I.ane  <h.st,ance.s  {D/t,  Df,  Dfi)  related  to  the  topo- 
ri.o  'iie  liiaji. 


Reproduced  From 
Best  Available  Copy 


2-7 


In  Figure  2.4,  the  transmitter  and  recdTer  are  bcated  atop  surface  heights  hj 
and  hn  respectively,  which  are  measured  from  a  sea  levd  xy-plane.  By  letting  the 
difference  in  height  be  Sha  =|  hr  —  h/i  |  ,  Rb  can  be  expressed  as: 


Rb  =  y(«hfl)>  +  (£>B)»)  (2.4) 

Similarly  Db  can  be  expressed  in  terms  of  Rb  and  Shs  as: 


DB  =  sl{RB?-(ShBf 

(2.5) 

By  defining  Shr  =|  hr  —  he  |  and  Shu  =|  h^  —  he  |  where  he  is  the  height  of  the 

dutter  at  the  point  of  reflection,  expressions  similar  to  Equations  2.4  and  2.5  can  be 

derived  for  Rt,  Rr,  Dj,  and  Dr. 

Rt  —  -f-  {Dry 

(2.6) 

Dt  =  yl(RT?  -  (^hr)’ 

(2.7) 

(2.8) 

DR=yJ[RRf^{8hRf 

(2.9) 

A  plot  of  the  equirange  contour  lines  is  illustrated  in  Figure  2.5.  The  concentric 
dlipses  have  the  transmitter  and  receiver  as  focal  points.  For  any  range  gate  interval, 
it  is  possible  for  the  target  dutter  cell  to  indude  any  of  the  area  between  the  cor¬ 
responding  equirange  dlipses.  The  goal  of  this  research  is  to  determine  the  bistatic 
dutter  power  from  the  topography  of  the  area  surrounding  the  system.  Therefore, 
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Figure  2.5.  eqnirange  oontoor  Enes  for  &  ground  l>ased  bistatic  system. 

a  prcoesi  is  needed  to  identify  any  or  all  of  the  eqnirange  oontonrs  for  a  given  radar 
qrstem  on  a  topographic  map  of  the  area  surrounding  the  system. 

Figure  2.6  was  partially  extracted  from  Beyer  (4:198)  and  illustrates  the  basic 
geometry  for  any  ellipse,  fbr  an  dlipse  centered  at  point  (h^)  with  its  major  axis 
parallel  to  the  z  axis  of  a  cartesian  coordinate  system,  Beyer  (4:199)  lists  the  following 
expression  for  the  ellipse: 


(x-h)»  (y-h)» 

a*  62 


(2.10) 


Expressing  a,  6,  h,  and  k  from  Equation  2.10  in  terms  of  (he  known  system 
parameters  tr,  Dr,  Dt,  Dr,  and  n  will  allow  the  ddineation  of  the  eqnirange  con¬ 
tour  Enes  on  the  topographic  map.  Ibr  convenience  we  dioose  a  coordinate  system 
centered  on  the  transmitter  with  the  z  axis  oollinear  with  the  slant  range  Rr. 

Ibr  any  elEpse,  the  sum  of  the  distances  from  any  point  on  the  ellipse  to  the 
foci  equals  the  distance  2a  in  Figure  2.6  (4:199).  Applying  this  principle  to  Ugore  2.7 
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vbere  the  distances  from  any  p<^t  on  the  eqniiange  oontonr  to  the  fixs  ate  Rt  and 
Rr  yields: 


Rt  +  Rr 
2 


(2.11) 


Erom  Equation  2.3  it  was  shown  that  Rr  +  Rr-V  eqnal  to  Rr  +ncTR,  Substi¬ 
tuting  this  into  Equation  2.11  yields: 


a 


Rb  ncTB 
2 


(2.12) 


Beyer  (4:199)  states  that  the  distance  from  the  center  cf  the  dlipse  to  either 
fed,  de/,  is  given  by: 


de/  «  va*  —  ^ 


(2.13) 
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Figure  2.7.  Geometry  for  the  case  inien  Rt  —  Rr. 

i 

Rom  Figure  2.6  it  it  dear  that  the  distance  dgf  it  equal  to  /2b/2.  Substituting 
Rb/2  into  Equation  2.13  for  <4/  *nd  the  right  hand  side  of  Equation  2.12  for  a, 
aimplifying  and  solving  for  6  yields:  | 


fncTR{2RB  +  ncTR) 


(2.14) 


Substituting  these  expressions  for  a  and  b  ^uen  by  Equations  2.12  and  2.14 
into  Equation  2.10  and  noting  that  h  =  0  for  the  geometry  shown  in  figure  2.6 


(Rb  +  ncTRy  ncTR{2RB  +  nerR) 
Rearranging  terms  in  Equation  2.15  and  solving  for  y  yields: 


nerR(2RB  +  ncrji)[(AB  +  ncrR^  -  4(*  -  h)^] 


2(Rb  +  ncTB) 


Using  Equation  2.16  the  magnitude  and  direction  of  Rr  can  be  determined  by: 


Rt  =  yjy^  + 

(2.17) 

y 

Onr  =  arctan  - 

X 

(2.18) 

Substituting  Equation  2.17  and  Shr  =|  hr  —  he  |  into  Equation  2.7  yields 
an  equation  which  can  be  used  to  determine  an  exact  value  £>r  Dt  which  in  turns 
identifies  two  points  on  an  equirange  contour  line  on  the  topographic  map. 

DT  =  y/x^  +  y^-ihT-he)^  (2.19) 

Using  Equation  2.19  to  identify  equirange  contour  fines  on  a  topographic  map 
can  become  a  very  time  consuming  process  because  use  of  the  equation  is  iterative  in 
nature.  Fbr  a  selected  value  of  x,  a  corresponding  y  is  calcnlated.  These  values  along 
with  the  known  transmitter  site  height,  hr,  can  then  be  used  with  Equation  2.19  to 
determine  Dt  from  the  dutter  heights,  he,  along  the  direction  ^ven  by  Equation  2.18. 
The  process  is  a  trial  and  error  procedure  when  completed  manually  in  search  of  the 
values  of  Dt  and  he  which  satisfy  the  equation.  If  digitized  surface  height  data  were 
available  the  process  could  be  automated,  but  this  is  not  the  case. 

An  alternative  is  to  establish  an  acceptable  kvd  of  error  in  Dt  and  identify 
the  conditions  where  Dt  can  be  approximated  by  Rt.  Rt  u  easily  identified  and 
can  be  plotted  directly  on  the  topographic  map  hj  using  Equations  2.16. 

Let  the  error  assodated  with  assuming  that  Dr  can  be  approximated  by  Rt 
be  denoted  by  tRo  where  crd  ia  defined  as: 


The  miuriinmn  ifidth  of  any  dntier  oell  viU  occur  abng  the  major  axis  of  the 
eqnirange  ellipses  which  define  the  odl.  Fbr  an  arbitrary  datter  cdl,  the  distance 
travded  by  the  wave  from  the  transmitter  to  the  leading  edge  of  the  dlipse  defining 
the  cdl  is  c(7o  +  (n  —  l)rA/2).  Similarly,  the  distance  tr&vded  the  wave  to  the 
ontermost  dlipse  which  defines  the  dutter  cdl  is  e(To  +  nr/t/2).  Therefore  the 
mtorimnin  width  of  the  dutter  cdl  is  pven  by: 

>y.  =  e(r.+^)-c^r.+  ^~te^  (2.21) 

The  criterion  for  determining  the  quality  of  the  Dt  Rt  assumption  will 
be  how  much  error  is  introduced  into  the  width  of  the  dutter  cdl  and  in  turn  the 
location  of  the  dutter  cell  ddimiting  ellipses.  Letting  the  P  represent  a  numerical 
percentage  expression  the  condition  for  the  assumption  is  stated  as: 


erd  <  PWc  (2.22) 

Substituting  Equations  2.21,  2.20,  and  2.7  into  Equation  2.22  yidds: 


Substituting  Skr  =  ^(Ar  ~  he)’  into  Equation  2.23  and  solving  for  he  yidds: 


he  <  cT/i^ P  [2/1  —  P{ctrY\  +  hr 


(2.24) 


By  spedfying  a  percentage  of  the  total  range  cdl  width,  P,  as  an  allowable  error. 
Equation  2.24  can  be  used  to  identify  a  range  of  dutter  hdghts,  Aei  which  allow  Dt 
to  be  approximated  by  Rj.  ibr  areas  which  meet  the  conditions  established  by 
Equation  2.24,  Equations  2.16  and  2.17  can  be  used  to  identify  the  perimeters  d 
Q>ecific  range  cdls. 
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Figure  2.8.  Gutter  Blocks  Created  by  the  Antenna  Patterns. 

2.S.2  The  Antenna  Patterns.  Once  the  equirange  contour  lines  have  been 
identified,  the  actual  shape  of  the  target  area  which  contributes  to  the  dutter  power 
detected  at  the  receiver  can  be  determined  by  overlaying  the  the  qrstem  antexma 
patterns  as  illustrated  in  figure  2.3.  In  this  research  project,  the  nuun  beam  pat¬ 
terns  as  well  as  the  significant  ade  lobe  patterns  will  be  consideted.  In  Hgure  2.8 
Ormi  and  0Rmi  represent  the  3  dB  beamwidths  of  the  transmitting  and  receiving  an¬ 
tennas  nuiil  lobes.  The  notation  used  is  similar  fer  the  side  lobes.  The  side  lobe 
beamwidths  represented  by  0t»i  and  0jui  are  determined  by  the  width  of  the  side 
lobes  at  points  where  thdr  power  has  decreased  to  3  dB  below  maximum.  The 
dashed  lines  represents  the  equirange  contour  lines  which  delimit  the  range  cell. 

The  impact  of  including  the  side  lobe  patterns  is  that  a  particular  range  cell 
may  include  up  to  nine  individual  dutter  blocks  which  act  as  sources  of  scattered 
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energy  as  opposed  to  a  single  dutter  block  when  only  the  main  beams  axe  considered. 
The  additional  dutter  blocks  are  a  result  of  the  interaction  of  the  main  and  side  lobes 
of  the  two  antennas.  E2ach  antenna  pattern  intersection  in  a  single  range  cell  will  be 
treated  as  a  separate  dutter  cell  as  shown  in  Figure  2.8.  Here  the  individual  dutter 
cdls  are  labeled  al,  a2,  a3,  and  <i4  to  identify  the  different  possible  combinations  of 
gains  and  beamwidths. 

Li  the  diagram,  al  represents  the  intersection  of  both  main  beam  patterns 
while  a2  identifies  the  intersection  of  the  transmitters  main  lobe  and  the  receivers 
dde  lobes.  Similarly,  a3  identifies  the  receiver  main  lobe  to  transmitter  side  bbe 
intersection  and  a4  the  side  lobe  to  side  lobe  intersections. 


The  number  of  additional  dutter  bbcks  to  be  oonsideted  will  depend  on  the 
width  of  the  range  cell,  the  bistatic  angle,  the  base  fine  distance  between  the  trans¬ 
mitter  and  receiver  and  the  individual  antenna  patterns.  Note  that  in  Figure  2.8  the 
range  cell  shown  has  only  eight  dutter  blocks.  This  situation  will  not  arise  fi>r  the 
system  used  in  this  research  because  of  the  narrow  transmitter  beam  width  and  the 
short  distances  between  the  i'ansmitter,  receiver,  and  target  area. 


Each  dutter  cell  may  have  to  be  fiixther  subdivided  into  blocks  and  analyzed 
individually.  The  analysis  of  the  target  area  terrain  involves  identifying  the  heights 
those  areas  which  exhibit  similar  roughness  characteristics  and  surface  cover.  Each 
dmilar  block  of  surface  area  was  analyzed  to  determine  the  mean  and  variance  of 


the  heights  and  a  height  distribution  whidi  best  describes  the  block.  Ftem  these 
diaracteristics  it  was  possible  to  calculate  the  RCS  for  each  target  area  bloddy  For 
those  dutter  cells  which  engulf  more  than  one  fype  of  surface  area,  the  cell  w^  be 
fiuther  subdivided  into  dutter  blocks  corresponding  to  the  respective  target  Wea 
blocks.  The  RCS  horn  each  target  area  block  can  then  be  used  with  the  area  of  the 
dutter  block  to  determine  the  power  contribution  from  that  dutter  block.  The  total 
power  from  the  dutter  cell  will  be  the  sum  of  the  powers  from  the  individual  dutter 


blocks. 
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The  notation  used  in  Figure  2.8  is  not  meant  to  imply  that  all  clutter  cells 
labeled  similarly  will  have  the  same  type  of  height  distribution  or  even  identical 
physical  areas.  This  requires  that  the  receiver  dutter  power  formula  given  by  Equa¬ 
tion  2.1  be  modified  to  replace  a  single  transmitter  gain,  receiver  gain,  and  bistatic 
cross  section  with  summations  as  shown  below  to  account  for  the  multiple  gain  and 
cross  section  products. 


{pTmGRm<fh,\  +  +  GjtGRm  I^j=l  +  GtsGrj  IZt=l  <^ba4k) 

(2.25) 


Where: 

Grm  =  gaiii  in  transmitter  main  lobe 
Gtm  =  gain  in  transmitter  side  lobe 
GRm  =  gain  in  receiver  main  lobe 
Gr$  =  gain  in  receiver  side  lobe 

=  CToss  section  of  at‘*  dutter  cell 


2.4  Determining  the  Incoherent  Contribution  Power. 

2.4.1  Rough  Surface  Models.  Meiny  models  exist  which  describe  scattering 
from  rough  surfaces  (3:70-99),  (13.-671-753),  (9:75-89).  The  most  general  of  these 
is  presented  by  Bartick  (13:703-729).  He  addresses  two  specific  cases  for  incoherent 
scattering,  slightly  rough  surfaces  and  very  rough  surfaces.  The  approach  varies 
for  determining  the  normalized  RCS  depending  upon  the  roughness  of  the  surface. 
Bartick  denotes  the  scattering  coefficient  as  7(^,-,0,,^,)  (13:672)  which  he  defines  as: 


A 


(2.26) 


Where: 


2-16 


A  =  The  area  of  the  terrain  surface  being  characterized. 

=  The  average  BCS  of  the  surface  in  question. 

It  is  apparent  from  Equation  2.26  that  Barrick’s  7(^1,  4>t)  ^  often  denoted  as 

(7°  the  normalized  Radar  cross  section  (NRCS)  other  authors.  For  a  ^ven  target 
area,  Ekjuation  2.26  will  have  to  be  used  £>r  each  section  of  the  target  area  that 
displays  different  height  or  sur&oe  cover  diaracteristics.  The  target  area  composed 
of  homogeneous  dxaracteristics  used  to  determine  the  RCS  of  the  area  may  not  be 
totally  illnminated  by  the  wave  front  £>r  a  single  range  cell.  The  RCS  in  Equation  2.1 
is  the  total  cross  section  of  a  single  range  oell.  Therefore  ai,  for  a  given  range  cell  is 
given  by: 

=  =  (2.27) 

Where: 

Ac  is  the  area  of  the  terrain  illuminated  by  the  wavefront  which  corresponds 
to  a  given  range  oeU. 

Two  assumptions  stated  by  Barrick  (13:672)  are:  DtDr  >  A,  the  product  of  the 
target  to  transmitter/receiver  distances  must  be  greater  than  the  area  of  the  terrain 
diaracteiization  surface  in  question,  and  4^  A,  the  operating  wavdength  must  be 
much  less  than  the  area  of  the  terrain  characterization  surface.  These  assumptions 
can  be  used  to  set  limits  on  the  size  ci  the  homogeneous  areas  of  terrain  used  to 
develop  the  terrain  statistics  (A^  '<  DtDr). 

The  toughness  scale  is  determined  ly  comparing  the  tm*  height  of  the  surface 
(h)  to  the  wavdength  (A)  of  the  agnal  being  transmitted.  For  cases  where  h  <C  X, 
the  surface  is  considered  to  be  slightly  rough.  When  h  >  A,  the  surface  is  considered 
as  a  very  rough  surface. 


S.4.2  Incoherent  Scattering  Fivm  a  Slightly  Rough  Surface.  Baiiick  (13:703- 
706)  explains  incoherent  scatteiin>j  from  a  slighdy  xongh  homogeneous  surface  as: 


Ipq  =  CCS-*  9i  cos*  $,  \oc„\^  I 


(2.28) 


Where  p  and  q  represent  the  polarization  state  of  (he  incident  and  scattered 
waves.  In  Equation  2.28  the  quantity  I  was  defined  as  finiows: 

^  =  jT rp(r)Jo  {rkoy/^+(^  dr  =  ^ir  {koy'Q  +  (i)  (2.29) 

hi  Equation  2.29  p(r),  and  r  are  defined  as; 


4c  =  sin  0,-  —  sb  cos 


(2.30) 


4y  s  an  sin 


(2.31) 


^  (C(».y)^y,y0> 


(2.32) 


r  =  y/{x  -  z')*  +  (y  -  y')* 


(2.33) 


Equation  2.32  is  the  surface  height  correlation  coefficient,  hi  order  to  use 
Equation  2.29  the  distribution  of  the  surface  heights  must  be  known  or  a  height 
distribution  must  be  assumed.  The  two  most  commonly  used  distributions  (13:704) 
are  Gaussian  and  two  sided  exponential.  A  Gaussian  surface  height  distribution 
will  be  assumed  for  the  work  accomplished  during  this  research.  Based  on  this 
assumption,  Equation  2.32  becomes: 
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Substituting  Equation  2.34  into  Equation  2.29  and  carrying  out  the  integration 
yields: 


I  =  ir/*exp 


4 


(2.35) 


Where  the  quantity  /  in  Equations  2.34  and  2.35  represents  the  correlation 
distance,  the  distance  required  to  cause  the  correlation  function  to  reach  a  value  of 
1/c. 


The  final  term  required  fisr  Equation  2.28  is  an  expression  for  the  a„  term.  As 
previously  stated,  only  linear  polarizations  are  to  be  considered  here,  therefore  the 
only  terms  which  apply  to  the  task  at  hand  are  aju  uid  0(w.  The  terms  are  defined 
by  Barrick  (13:706)  as: 


Othk  = 


{Hr  - 1)  sin  Bi  daBt-oos  -  an*  Biyjuih  -  am*  B^  +  (^  ”  1) 

cos  Bi  +  ^Crlir-dn^B,^  008  B,  +  yj Crtlr  -  dv}  B^ 

(2.36) 


otw,  = 


(Cr  -  1)  an  Bj  sin  g,  -  cos  -  sin^  Bjyjcrlh  —  B^  -f  ^  {nr  —  1)  cos  <f>. 


^C,  cos  Bi  +  yJcrltr—dsi^Bi^  I^Cr  COS  B,  +  yj Crlh  -  ds}  B^ 


(2.37) 


Where: 


Cr  =The  relative  permittivity  of  the  target  area. 
Hr  =The  relative  permeability  of  the  target  area. 


Equations  2.36  and  2.37  can  be  amplified  by  assuming  that  objects  in  the 
target  area  will  probably  consist  of  nonmagnetic  materials.  This  implies  that  /^r  ^  1 
and  Equations  2.36  and  2..37  can  be  recast  as: 


Cthh  = 


(gf  -  1)  cos  <l>t 


^008  6i  +  —  sin^^  ^oos  0,  +  —  sin*^ 


(2.38) 


Q!»ti  = 


(Cr  —  1)  sin 6ian0t  —  cos ^,^c,  —  sin^ 0i\/^  —  sin^ 6,J 
^c,  cos  0i  +  \/er  —  an^^ij  cos  —  an*^ 


(2.39) 


The  incoherent  scattering  from  slightly  rough  surfaces  typically  is  caused  by 
the  material  covering  the  terrain  surface  like  grass,  crops,  and  foliage.  A  detailed 
description  of  these  items  is  required  to  develop  a  statistical  characterization  of  their 
heights  and  composition.  Unfortunately,  this  type  of  information  is  not  available 
from  topographic  maps. 

The  maps  do  provide  a  general  description  of  the  type  inland  cover  in  categories 
Eke  forest,  swamp,  sand  etc.  More  spediic  information  about  the  ground  cover  is 
required  to  implement  the  procedure  outlined  in  this  section  therefore,  incoherent 
scattering  from  sEghtly  tough  surfaces  will  not  be  demonstrated  in  this  report. 

2.4.S  Incoherent  Scattering  From  Very  Rough  Surfaces  As  the  roughness  of 
the  surface  increases,  the  scattered  power  attains  more  contribution  from  the  in¬ 
coherently  scattered  power  and  less  from  the  coherently  scattered  power.  Barrick 
states  (13:719)  for  koh  >  5  the  scattered  power  is  essentially  incoherent. 

The  approach  for  deriving  the  very  tough  surface  solutions  differs  from  the 
slightly  rough  surface  iq)proach.  The  principle  of  physical  optics  was  used  by  Bar¬ 
rick  to  solve  the  Stratton-Chu  integral  equation  to  describe  the  scattering  from  very 
rough  surfaces.  His  derivation  and  analysis  are  based  on  the  tangent  plane  approxi- 
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•  mation  to  tl;r  stufaco  otrrrnt  prrKcnt  on  n  iluite  (»octinn  of  torrain  to  .solve 

the  integral  c<]nation.  A  stalisliral  roprrsoiit.ation  of  the  randomly  oiiontc'd  t'’rr.'un 
stirface  ludglits  is  rcqnirrj  to  ojcprcss  the  scattering  .as  simcular  reilection  frcjin  the 
raisdom  surfarc'j?  in  the  targ.et  jura  (13:711}).  llis  approach  (13:720)  invokes  the 
following  restrictions; 

1.  The  radius  of  curvature  A. 

2.  The  rough  ne.ss  is  Isotropic. 

3.  The  correlation  length  / 

4.  The  (dTects  of  multij»lc  scattering  .and  shallowing  were  neglected. 

Ihosed  upon  these  a.ssuniptions,  the  average  NliCS  .as  derived  by  Barrick  (13:720) 
is  given  by: 

(2-10) 

In  Erpiation  2,10  and  J  .are  .uiaio  uiiis  to  the  quantities  Op,  and  7  for  the 
slightly  rough  Mirfaee.  Whc'rc  J  is  defined  by  Barrirk  .as; 

J  =  'X  JJ'  rJ„  (rho^a  +  (j)  (i  Ao^.,  -iUr,r)  dr  (2.11) 

\Vliere; 

7o(i)  i.s  a  ;vro  order  cylindrical  Ik  s  cl  function. 

^  and  (^'  are  .siirface  height  r.andoni  v:iri,ables  .'oparated  a  di.st.ancc  r. 

Is  t!ie  joint  characteristic  function  of  the  .surface  heii'ht  random 

variables. 

Tl*i5  joint  cJiaracteri.stic  function  i.s  the  I'luirier  transfonu  of  the  joint  prob.a- 
bility  density  function  (.IPDF),  /  j(C,C’).  v.-ldch  de.scribes  the  surface  heights.  As 
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with  the  I  function  from  Equation  2.29  for  the  ilightly  rough  surface,  the  density 
function  describing  the  heights  must  be  Imown  or  assumed  to  calculate  a  value  ioz 
the  J  function.  This  JPDF  will  again  be  assumed  to  be  Gaussian  and  given  by: 


mx) 


2ff  h*  (1  -  p*(r))  ^  ^  I  2A»  (1  -  J 


(2.42) 


la  Equation  2.42,  (  and  represent  the  surface  heights  at  two  locations  (z,  y) 
and  (x',y')  and  p(r)  represents  the  Gaussian  correlation  fimction  ^ven  by: 

p(r)  =  exp  (2.43) 

Where: 

r*  =  (z  —  z')’  +  (y  —  yO*  Th®  distance  between  two  points. 

P  =  The  square  of  the  correlation  distance. 

Barrick  (13:721)  argues  that  the  Gaussian  correlation  coefficient,  Equation  2.43, 
can  be  expressed  the  first  two  terms  of  a  power  series  expansion  for  his  develop¬ 
ment  of  the  high  frequency  scattering  modd  as: 


p,  » 1  - 


(2.44) 


Using  the  Gaussian  suri[ace  height  JPDF  ^ven  in  Equation  2.42  and  the  small 
signal  ^>proximation  fi3r  the '^Gaussian  surface  height  oorrdation  function  given  in 
Equation  2.43  in  Elquation  2.4^  yields  the  following  expression  after  integration: 

Where  the  quantities  in  Equation  2.45  s’,  (y,  and  ^  have  been  defined  as: 


2-22 


=:  on  —  tin  oot 


(2.46) 


(y  =  an  $i  ain 


(2.47) 


(,  =  —  cog  tf,-  —  008 


(2.48) 


a  4A» 
^  ~  P 


(2.49) 


Where: 

is  the  mean  squared  roughness  height. 

/  it  the  surface  correlation  length. 

The  scattering  matrix  in  Equation  2.40  is  represented  by  where  again  the 
subscripts  p  and  q  denote  the  polarization  state  of  the  transmit  and  leoeiTe  antennas 
respectively.  In  this  research,  only  vertical  and  horizontal  linear  polarizations  were 
considered  and  the  scattering  matrix  dements  were  defined  as: 


P  _  oaa3J?ll(t)  +  sin  g,-  an  (f,  sin^  ^»ili.(0 

**  C1O4 


(2.50) 


_  -  an  9i  an  9,  sm^  ^>iZn(*)  - 


(2.51) 


The  expressions  iZx(*)  >>^<1  -^K*)  ^  modified  Resnel  reflection  coefficients 
given  by  the  following  expressions: 


e,  cos  t  -  y/€rllr  —  gp  » 
e,  cos  t  + 


(2.52) 
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The  angle  i  used  in  Equations  2.52  and  2.53  has  been  defined  as  follows  by 
Barrick: 

oos  t  =  —  an  BiwaBt  cot  + CO*  6i  cos  9,  (2.54) 

The  quantities  ui,  uj,  03,  and  04  used  in  Equations  2.50  and  2.51  have  been 
defined  as  follows: 


ui  =  1  +  sin  9i  sin  oos  —  oos  0i  cot  0,  (2.55) 

U]  =  008  an  +  >in  0i  cos  0,  oos  (2.56) 

03  =  sin  0i  an  9,  +  oos  sin  cos  (2.57) 

04  =  008  +  008  01  (2.58) 


By  using  the  expressions  ^ven  above,  the  NRCS  of  a  selected  block  of  ter¬ 
rain  which  meets  all  the  restrictions  fisted  can  be  calculated  £>r  any  given  bistatic 
geometry.  The  NRCS  can  then  be  multiplied  by  the  area  of  this  terrain  which  is 
illuminated  by  both  transmitting  and  receiving  antennas  fi>r  a  particular  range  cell 
to  determine  the  the  RCS.  Once  the  RCS  is  known,  the  scattered  power  from  the 
section  of  terrain  can  be  determined  by  using  Equation  2.25. 


III.  Application  of  the  Process  and  Results 

S.l  Known  System  Parameters. 

la  this  diapter  the  process  described  in  Chapter  11  vrill  be  ^>plied  to  a  ground 
based  Ustatic  qrstem  operated  by  the  sponsor,  RADC/EECE.  The  process  dements 
which  can  be  implemented  from  information  obtained  from  topographic  maps  will 
be  illustrated  and  results  presented.  The  starting  point  is  to  identify  all  known 
parameters  for  the  qrstem  in  question. 

The  transmitting  hardware  is  a  pulsed  3  GHz  S  band  weather  radar  ^stem  lo¬ 
cated  near  Stow  Massachusetts  and  operated  by  the  Air  Esrce  Geophysics  Laboratory 
(AFGL).  The  transmitter  radiates  a  peak  power  of  800  KW  using  a  28  foot  parabolic 
dish  antenna.  The  antenna  is  located  atop  a  15  foot  tower  on  a  hill  top  which  has 
an  devation  of  ^>proximatdy  93  meters,  therefore  the  quantity  hr  ^  97.572  meters. 
The  transmitter  has  a  PRP  of  1  KHz  and  a  pulse  width  of  1  ftS.  The  beam  width  of 
the  antenna’s  main  lobe  was  estimated  by  the  qmnsor  to  be  approximately  1*. 

The  leceiTing  hardware  is  bcated  atop  Prospect  £511  in  Waltham  Massachusetts 
and  employs  a  standard  gain  20  dB  horn  antenna.  Prospect  ffill  is  approximatdy 
144  meters  above  sea  levd,  therefore  the  quantity  kn  fts  144  meters.  The  receiver  is 
capable  of  operating  with  a  variety  of  range  gate  intervals  but  will  use  a  1  /xS  interval 
iat  this  project,  therefore  the  quantity  rn  =  IftS. 

The  linear  distance  between  the  transmitter  and  receiver  sites  is  19.248  Km 
and  the  baseline  slant  range  am  be  determined  from  Equation  !2.4  using  hr  =  97.572, 
hfi  =  144,  and  Db  =  19248.  Performing  the  calculation  yields  Rb  Db  =  19248 


meters. 


S.2  Establishing  the  Clutter  Cells. 

Ib  establish  the  clutter  odls  three  patterns  most  be  identified.  The  equirange 
contour  Hues  must  be  identified  on  the  topographic  map.  hi  addition  both  antenna 
patterns  must  be  identified.  The  intersection  of  the  antenna  patterns  and  a  particular 
range  cell  identifies  the  clutter  cell  dimension  ^plicable  to  that  particular  range  cell. 

3.2.1  Identifying  the  Equirange  Contour  Lines.  The  first  step  is  to  determine 
the  equirange  contour  Hues.  Both  sites  are  located  on  a  sin^e  USGS  topographic 
map,  The  Maynard  Massachusetts  7.5  by  15  minute  map.  A  rectangular  coordinate 
system  uras  established  on  the  topographic  map  such  that  the  transmitter  site  was 
located  at  the  origin  and  the  base  line  between  the  transmitter  and  receiver  was 
taken  as  the  x  axis.  The  coordinate  system  was  scaled  so  that  1  cm  corresponds  to  an 
actual  distance  of  250  meters.  A  fortran  program  was  written  based  on  Equation  2.10 
which  produces  y  values  for  250  meter  increments  of  the  x  coordinate.  The  equirange 
contour  lines  were  manually  plotted  on  the  map  for  the  first  22  range  cells. 

The  limitations  initially  placed  on  the  data  collection  area  by  the  sponsor 
were  that  the  powers  from  3  to  5  range  cells  were  to  be  estimated  with  the  bistatic 
angle,  held  between  120”  and  160”.  These  requirements  were  altered  during  the 
course  of  the  research  in  that  the  sponsor  specified  three  actual  target  sites  on  the 
aforementioned  map. 

3.2.2  Identifying  the  Antenna  Patterns.  little  information  was  available  about 
the  antenna  diaracteristics  in  addition  to  the  parabolic  dish  diameter  and  the  gain 
of  the  horn  antenna.  However  from  the  information  available  it  was  necessary  to  de¬ 
termine  the  beam  widths  of  the  main  lobes  and  side  lobes  of  both  antennas  and  their 
locations  relative  to  the  antenna  bore  sight  so  that  the  dutter  block  area  determined 

the  intersection  of  the  two  patterns  could  be  determined. 
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.  S.2.2.1  The  Thinsmitting  Antenna.  As  previously  stated  the  only  in¬ 
formation  available  on  the  transmitting  antenna  pattern  vvas  that  the  main  lobe 
beam  tvidth  'was  approximatdy  1".  In  order  to  identify  the  side  lobe  beam  width 
and  relative  locations  of  the  side  lobes  as  wdl  as  the  gains  associated  with  the  beam 
widths  it  was  necessary  to  perform  some  reverse  engineering.  No  information  was 
available  to  the  researcher  about  the  dish  feed  antenna  or  the  illumination  pattern. 

As  a  first  attempt  to  identify  a  feed  pattern  which  produced  a  1”  beam-width, 
the  parabolic  reflector  was  treated  as  uniformly  illuminated  circular  aperture,  where 
the  transmitting  antenna  had  a  diameter  of  8.5344  meters.  Kraus  (7:344)  gives  the 
normalized  Add  pattern  as: 


£.(«)  =  2 


Jx 


(’f)me 


Where: 

3s  the  diameter  of  the  aperture 
$  as  the  angle  with  respect  to  the  normal  to  the  aperture 
Ji  first  order  Bessel  fonction 

The  normalized  power  pattern  can  be  obtained  from  Equation  3.1  as: 


(3.1) 


The  beam  width  for  the  main  beams  -will  be  taken  as  the  3  dB  beam  width  of 


the  normalized  antenna  power  patterns.  The  3  dB  points  occur  at  the  point  where 
Equation  3.2  drops  to  a  value  of  1/2  from  the  mATimnni  power  point  which  occurs 
on  bore  sight.  Letting  x  s=  wZl/Asin  9,  implies  that: 


1 

X  2y/2 


»  0.35355 


(3.3) 


Uuing  Balanic’  Ji(a:)/x  table  (1:943)  and  interpolating  to  find  the  value  of  x 
for  Ji{x)fx  =  0.35355  yields  x  =  1.61612.  But  x  =  irDfXdnO,  therefore  the  3  dB 
point  occurs  at  0  =  arc3iu(1.61612A/jri)).  The  half  power  beam  width  (HPBW)  is 
equal  to  twice  the  3  dB  point,  hence: 

«  .  /1.61612A\ 

HPBVV^t  =  2  arcsin  (  — — — j  (3.4) 

Applying  Equation  3.4  to  the  transmitting  antenna  yields  a  beam  width  of  0.7®. 
In  an  attempt  to  locate  a  feed  pattern  which  produced  more  acceptable  results,  the 
feod  rattern  was  assumed  to  have  a  parabolic  taper.  Stutzman  and  Thiele  (17:420- 
‘iil)  provide  information  in  tabular  format  for  several  orders  of  parabolic  taper 
and  edge  illumination.  Using  the  HPBW  as  a  guide  it  was  determined  that  the 
transmitter  pattern  more  closely  matches  a  second  order  parabolic  taper.  They 
estimate  the  HPBW  in  radians  as: 

HPBlV{rad)  =  —  =  —  (3.5) 

Using  a  wavelength  of  10  cm  and  a  dish  diameter  of  853.44  cm  in  Equation  3.5 
yields  a  HPBW  which  is  approximately  1®. 

HPBW{rad)  =  =  0.0172244rad  =  0.986®  (3.6) 

The  normalized  field  pattern  associated  with  a  second  order  parabolic  taper  is 
given  by: 


En{0)  = 


sin  0) 

(^sin0)3 


(3.7) 


The  normalized  power  pattern  of  the  antenna  is  equal  to  the  square  of  the 
normalized  field  pattern,  therefore: 
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The  location  of  the  side  bbes  and  the  lespective  3  dB  down  beam  widths  must 
be  determined  £br  the  antenna  pattern.  The  first  side  bbes  occur  at  the  first  peak 
of  the  Bessel  function  away  from  bore  sight.  Beferring  the  the  Balauis  tables  again 
(1:942)  this  peak  at  x  »  7.9167.  But  x  =  rDI\ma.0,  therefore  setting  the  two 
quantities  equal  and  rearranging  to  solve  br  0  yields  the  center  line  direction  br  the 
first  side  bbes. 


.  ^  .  /7.9167A\ 

6^  =  ±  axcsin  I  J  (3.9) 

Solving  Equation  3.9  to  determine  the  angle  which  corresponds  to  the  side  bbe 
center  lines  br  the  transmitter  (0tit)  yields  =  1>70”.  The  normalized  power  bvels 
at  the  side  bbe  center  lines  was  obtained  from  Equation  3.8  by  setting  x  =  7.9167 
and  solving  which  yielded  a  power  of  816.956  /tW  or  —30.9  dB.  At  the  3  dB  down 
points  from  the  side  bbe  center  line  this  power  is  cut  in  half  to  408.48  filV.  Therefore 
solving  br  the  corresponding  an^e  Fn(^)  =  408.48  fiW  yields  a  recursion  equation 
in  z.  Solving  the  equation  yidds  x  »  8.81.  Replacing  7.9167  in  Equation  3.9  with 
8.81  and  solving  identifies  the  bcation  of  one  of  the  3  dB  down  angles. 

=  (3.10) 

Using  Equation  3.10  to  find  the  side  bbe  3  dB  down  angles  br  the  transmitter 
produces  1.88”.  The  ade  bbe  HPBW  equals  twice  the  difference  between  the  side 
bbe  center  line  angles  and  the  3  dB  down  angles.  &nce  br  the  transmitter: 


HPBW^  =  2(0,u  -  W)  =  0.37” 


(3.11) 


The  final  process  in  determining  the  transmitter  antenna  pattern  was  to  de¬ 
termine  the  appropriate  gains  fer  the  main  and  side  lobe  beam  widths.  Barksdale 
(2:1-54)  identifies  the  directive  gain  fiinction  of  an  antenna  as: 

D,i$)  =  GoFnie)  (3.12) 

Where  Go  in  Eiquation  3.12  is  the  ma-riTnnm  power  gain  of  the  antenna.  Equa¬ 
tion  3.12  has  been  modified  to  no  ^  variations  since  the  parabolic  antenna  in  question 
has  antenna  patterns  which  are  qrmmettic  in  the  ^  direction.  The  maximum  power 
gain  of  the  antenna,  Go,  is  equal  to  the  cadency  of  the  ant'^nna  times  the  directivity 
of  the  antenna. 


Go  =  cDo  (3.13) 

The  effidency  of  a  parabolic  reflector  (lOfl)  can  be  estimated  as,  e  =  0.55. 
Kraus  (7:345)  identifies  the  directivity  of  an  antenna  as: 


^  (3.14) 

Where  Ap  in  Equation  3.14  is  the  physical  area  of  the  aperture.  For  the 
parabolic  reflectors  the  aperture  area  is  the  area  of  a  drde  of  diameter  D.  Sub¬ 
stituting  Ap  =  s’(^)^  into  Equation  3.14  results  in: 


Do 


(3.15) 


Combining  Equation  3.15  and  the  effidency  estimate  of  0.55  in  Equation  3.13 
results  in  a  maximum  power  gain  of: 
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Applying  Equation  3.16  to  the  transmitting  antenna  in  thic  bistatic  qrstem 
yields  a  nm-TimiiTn  gain  for  the  transmitter  of: 

Go,  =  (0.55)  =  39537.4  =  45.97  dB  (3.17) 

Applying  the  results  d  Equation  3.17  to  Equation  3.12  results  in  gain  functions 
which  can  be  used  to  identify  the  gains  associated  with  the  main  and  side  lobe  besun 
widths  of  the  antenna.  The  gain  of  the  antenna  main  beam  was  taken  to  be  equal  to 
the  majrimuTn  directive  gain.  The  gain  for  the  side  bbe  beam  widths  were  calculated 
as  the  gain  at  the  side  lobe  center  angles,  0  =  9m  =  1.70°. 

DM  =  Go,FM  34.43  =  15.37  dB  (3.18) 

S.2.2.2  The  Receiving  Antenna  As  previously  mentioned  the  only  in¬ 
formation  provided  was  that  the  receiving  antenna  was  a  20  dB  standard  gain  horn 
antennsL  Ikom  the  literature  reviewed  the  term  standard  gain  was  most  often  ap¬ 
plied  to  pyramidal  horn  antennas  designed  for  optimum  gain.  Though  this  was  not 
explicitly  stated  in  the  material,  the  following  procedures  are  based  on  this  assump¬ 
tion. 

Stutzman  and  Thiele  (17:411-415)  outline  a  procedure  for  designing  an  opti¬ 
mum  gain  pyramidal  horn.  The  results  of  their  design  procedure  are  the  pertinent 
physical  dimensions  of  the  antenna.  Tbese  dimensions  are  illustrated  in  figure  3.1. 
fVom  these  dimensions  the  gains  and  beam  widths  can  be  estimated  from  the  uni¬ 
versal  radiation  patterns  provided  (17:409). 

Ikom  the  conditions  for  optimum  gain  of  E-plane  and  H-plane  sectoral  horns, 
Stutzman  and  Thiele  (17:413)  develop  a  recursion  formula  for  A»itigning  an  optimum 
gain  pyramidal  horn  illustrated  in  Equation  3.19. 
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Figure  3.1.  Fjrramidal  Horn  Geometry. 
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(3.19) 


(2<t-1) 


.  185r^<r 


Equation  3.19  can  be  used  to  solve  for  <r  provided  the  waveguide  dimensions  a 
and  b  feeding  the  horn  and  the  gain  G  are  known.  In  Equation  3.19  cr  is  gi  —  jy: 


<7  — 


A 


(3.20) 


Therefore  Equation  3.19  can  be  used  to  solve  for  a  which  in  turn  can  be  used 
to  determine  the  dimension  i£.  With  the  value  of  Ie  Equation  3.21  can  be  used  to 
determine  B. 


B  =  ^J2\R2  «  (3.21) 

With  a  value  for  B  Equation  3.22  can  be  used  to  determine  A. 

G  =  (3.22) 

Where  tap  in  Equation  3.22  is  the  aperture  efficiency  of  the  anteima.  An 
aperture  efficiency  of  0.5  is  a  common  value  used  for  horn  antennas.  The  quantity 
Ap  in  Equation  3.22  represents  the  physicrJ  area  of  the  antenna  aperture  which  is 
equal  to  the  product  AB.  With  the  lengths  A  and  B  established,  the  universal 
radiation  patterns  provided  by  Stutzman  and  Thiele  (17:409)  can  then  be  used  to 
identify  the  beam  widths  of  the  main  and  side  lobes  for  both  the  E-plane  and  H-plane 
sectors  as  well  as  the  respective  directions  of  the  side  lobes  relative  to  bore  site. 

The  waveguide  dimensions  were  obtained  from  Pozar  (11:716),  where  he  iden¬ 
tifies  a  WR-284  waveguide  as  an  appropriate  choice  for  S  band.  The  dimensions 
for  WR-284  are  a  =  7.214  cm  and  b  =  3.404  cm.  Stutzman  and  Thiele  suggest 
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an  initial  value  of  (7  =  G/2Tr-^.  After  many  iterations,  it  was  determined  that 
<r  »  6.21430022.  Using  this  value  in  Equation  3.20  and  solving  for  t  yidds: 

i  =  \a  =  (.1)(6.21430022)  «  62.413cm  (3.23) 

Using  this  restilt  in  Equation  3.21  to  solve  fot  B  yidds: 


B  «  =  y(2)(10)(62.413)  =  35.33cm  (3.24) 

Using  this  value  &r  £  in  Equation  3.22  and  solving  &r  A  yidds: 


GX*  (100)(10)^ 

2vB  “  (2)(35.33)t 


45.0dcm 


(3.25) 


Using  these  values  and  the  universal  radiation  patterns  (17:404,409)  the  neces¬ 
sary  patterns  were  identified.  In  the  Erplane,  the  3  dB  beam  width  was  established 
as  14.14”.  The  first  side  lobe  occurs  at  23.6”  has  a  beam  width  of  7.6'^  and  is  9.1  dB 
dovni  from  the  main  lobe,  fri  the  E-plane,  the  3  dB  beam  width  was  established  as 
11.1”  and  the  side  lobes  are  not  distinguishable  from  the  main  beam. 


S.2.S  Identifying  the  Clutter  Cells.  The  clutter  cells  assume  the  shape  out¬ 
lined  by  the  liteisection  of  a  sdected  range  cells,  the  transmitter  antenna  patterns, 
and  lecdver  antenna  patterns.  Before  the  antenna  patterns  from  Section  3.2.2  can 
be  overlaid  on  the  equirange  contour  lines  which  delineate  a  ^edfic  range  cell,  the 
bistatic  angle  for  the  system  must  be  sdected.  As  previously  mentioned,  the  require¬ 
ments  were  dtered  by  the  sponsor  during  the  course  of  the  research.  Originally  three 
bistatic  angles  were  to  be  sdected  between  120”  and  160”.  In  light  of  the  fact  that 
an  adequate  shadowing  function  had  not  been  devdoped  and  use  of  a  shadowing 
function  would  require  digitalization  of  all  terrain  hdghts  between  the  transmitter 
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•ad  target  as  well  as  between  the  receiver  and  target,  three  sites  were  jointly  selected 
by  the  researcher  and  the  sponsor. 

Site  one  was  a  hiUy  forested  area  located  approximately  500  meters  south  of 
Lincoln,  Massachusetts.  The  distance  horn  site  one  to  the  transmitter  {Dj)  was 
measured  to  be  14.4625  Km.  The  distance  from  site  one  to  the  receiver  {Dr)  was 
measured  to  be  5.775  Km.  The  bistatic  angle  was  determined  by  applying  the  law 
of  cosines  to  the  geometry  involved,  'rb?  general  form  is  provided  in  Equation  3.25. 
Applying  the  values  for  Pr,  Dr,  and  the  basebne  distance  between  the  transmitter 
and  receiver  (Dr  —  19.248A'm)  to  Equation  3.26  yielded  a  histatic  angle  of  140.0° 
for  site  one. 


=  arccos 


2DrDt 


(3.26) 


Site  two  was  fhirhaven  Hill  located  approximately  2.5  Km  south  of  Concord 
Massachusetts.  The  distance  from  site  two  to  the  transmitter  (Pr)  measured  to 
be  11.050  Km.  The  distance  from  site  two  to  the  receiver  {Dr)  was  measured  to  be 
9.875  Km.  Using  these  values  in  Equation  3.26  resulted  in  a  bistatic  angle  of  133.7° 
for  the  physical  layout  of  site  two. 

Site  three  was  Pine  Hill  located  approximately  500  meters  west  of  the  Hanscom 
Air  Force  Base  runways.  The  distance  from  site  one  to  the  transmitter  (Pr)  was 
measured  to  be  16.1125  Km.  The  distance  from  site  one  to  the  receiver  {Dr)  was 
measured  to  be  9.9625  Km.  Substituting  these  values  into  Equation  3.26  yielded  a 
bistatic  angle  for  site  three  of  approximately  92.1”,  which  is  outside  of  the  original 
target  range  for  bistatic  angles. 

The  antenna  patterns  were  then  manually  plotted  on  the  topographic  map 
at  each  of  the  selected  ates.  The  intersection  of  the  transmit  and  receive  antenna 
patterns  was  used  to  identify  the  portions  of  range  cells  illuminated  at  each  site. 
This  intersection  of  the  transmit  and  receive  antenna  patterns  along  with  the  range 
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cell  constitutes  a  dutter  cell.  At  site  one,  the  transmit  antenna’s  main  lobe  and  side 
lobe  patterns  intersect  with  the  recdve  antenna’s  main  lobe  pattern  over  range  cells 
3,  4,  5,  and  6.  The  recdver  antenna  ode  lobe  pattern  intersects  with  the  transmit 
-antenna’s  main  and  side  lobe  patterns  over  range  cells  6,  7,  8,  and  9. 

Fbr  site  two,  the  transmit  main  and  side  lobes  intersect  with  the  receive  main 
lobe  over  range  cells  4,  5,  6,  7,  and  8.  The  receive  ante^a  side  bbe  intersects 
the  transmit  antenna’s  main  and  side  lobes  over  range  cdls  10,  11,  12,  13,  14,  and 
15.  As  the  bistatic  angle  decreases,  the  antenna  pattern  intersection  points  occur 
at  distances  farther  from  the  baseline  between  the  transmitter  and  receiver.  As 
the  intersections  occtu  hirther  from  the  baseline,  the  range  cells  get  progressively 
narrower  and  hence  more  range  cells  are  illuminated. 

Ear  site  three,  the  transmit  main  and  side  lobes  intersect  with  the  receive  main 
lobe  over  range  cells  18,  19,  20,  21,  22,  23,  24,  25,  26, 27,  28,  and  29.  The  receive  side 
lobe  intersects  with  the  transmit  main  and  side  lobes  over  many  range  cells  starting 
from  cell  number  24  and  continuing  off  the  map. 

Li  order  to  apply  Barrick’s  rough  surface  scattering  formulas,  it  was  necessary 
to  diattafterize  the  terrain  being  illuminated  in  the  dutter  cells  statistically.  The 
terrain  to  be  characterized  was  first  divided  into  rectangular  blocks  as  shown  in 
Figure  3.3  such  that  each  Hock  covers  an  area  of  terrain  with  homogeneous  diarac- 
teristics.  Unfortimatdy,  nature  rarely  simplifies  the  procedure  by  allowing  an  entire 
area  illuminated  to  be  composed  of  a  single  area  cf  homogeneous  characteristics.  An 
example  at  this  point  may  be  helpful. 

Assume  that  a  particular  dutter  cell  has  been  identified  on  the  topographic 
map.  The  terrain  area  of  the  cell  is  partially  composed  of  forested  hiUs  and  partially 
composed  of  tolling  grass  land.  The  site  would  have  to  be  divided  into  two  rect¬ 
angular  blocks,  one  to  characterize  the  forested  area  and  a  second  to  diaracterize 
the  grass  land.  The  intersection  of  the  dutter  cell  and  the  forested  characterization 
block  would  be  considered  a  dutter  block.  Similarly,  the  intersection  of  the  grass  land 


diaracterization  block  and  the  dutter  cell  ivould  comprise  a  second  dutter  block.  To 
find  the  scattered  power  from  this  section  of  the  range  odl  would  require  solving  frr 
the  dutter  power  from  each  dutter  block  individually.  The  total  power  from  the  two 
dutter  blocks  would  be  the  sum  of  the  individual  magnitudes. 

S.S  Terrain  Site  Characterization. 

It  is  easy  to  see  how  performing  these  calculations  manually  could  quickly 
become  an  accounting  nightmare.  The  area  illuminated  by  the  receive  antenna  main 
lobes  and  the  transmitter  antenna’s  total  pattern  covers  portions  of  five  range  cells. 
The  receive  antenna’s  side  lobes  intersect  with  the  transmitter  anteima’s  full  pattern 
over  portions  of  six  range  cells.  |  In  this  situation,  each  range  cell  is  composed  of  from 

j 

one  to  three  dutter  cells.  The  complexity  increases  with  the  number  of  dutter  blocks 
involved  with  each  dutter  cell.  I  Fbr  these  reasons,  only  three  range  cells  at  each  site 
were  sdected  for  the  demonstration  of  the  process  described  in  Chapter  two. 

i 

I 

Several  pieces  of  statistical  information  from  each  homogeneous  terrain  area 
to  be  diaracterized  were  required  to  iq>ply  Barrick’s  formulas  for  determining  the 
NRCS  of  the  areas  in  questionJ 

hirst,  recall  from  Section  2.4.3  that  knowledge  of  the  height  distribution  was 

required  to  solve  the  surface  coirrent  integral  in  Equation  2.41.  La  this  research,  the 

! 

distribution  has  been  auumed  to  be  a  joint  Gaussian  probability  density  fimction 
denoted  by  the  random  variable  Y).  One  advantage  to  this  assumption  is  that 
Barrick  (13:721)  has  provided  a  dosed  form  solution  for  Equation  2.41  when  the 
distribution  is  Gaussian. 

Barrick  also  provides  a  dosed  form  solution  for  a  joint  exponential  height  dis¬ 
tribution  (13:721)  however,  Roseman  (12:1-2)  disputes  the  validity  of  the  joint  ex¬ 
ponential  distribution  as  a  valid  joint  probability  distribution  fimction  and  raises 
questions  as  to  its  validity  for  representing  rough  surfaces.  Therefore,  for  any  distri- 
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7. 


Figoze  3.2.  Terrain  Scattering  Geometry. 

bation  other  than  the  joint  Ganssian  pdf,  Equation  2.41  would  probably  have  to  be 
aolved  numerically. 

Second,  the  reference  |dane  depicted  in  Figure  3.2  is  oo-planer  with  a  plane 
passing  through  the  mean  value  of  the  surface  height  distribution.  Therefore,  the 
sample  mean  of  the  random  variable  ({Xf  Y)  was  required  to  determine  the  incident 
angle  and  scattering  angles  of  the  electromagnetic  energy  at  the  terrain  site. 

The  third  statistical  quantity  required  from  the  surface  height  distribution 
was  the  sample  variance  of  the  surface  heights.  The  quantity  was  used  directly 
m  Barrick’s  dosed  form  solution  of  Equation  2.41  (13:721).  Baizick  refers  to  the 
variance  of  the  sample  mean  as  the  mean  square  roughness  height  and  identifies  it 
with  the  qrmbol  A^. 

The  final  data,  and  probably  the  most  difScult  information  to  estimate,  is  the 
correlation  length  of  the  joint  Gaussian  distribution. 

l^th  the  knowledge  of  the  these  four  items,  Barrick’s  formulation  can  be  used 
to  determine  a  value  of  NRCS  for  the  terrain  area  covered  by  the  statistical  distri- 


/ 
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bution.  Once  the  NRCS  was  known,  the  RCS  of  that  portion  of  the  terrain  visible 
simultaneously  to  both  antennas  was  easily  calculated. 


The  terrain  areas  to  be  characterized  were  divided  into  rectangular  areas  with 
homogeneous  diaracteristics  as  shown  in  Figure  3.3.  The  total  area  was  divided 
along  each  axis  into  n  blocks  of  width  Sf.  The  result  was  blocks  each  with  an 
area  of  6^  square  meters. 

The  terrain  height  at  any  point  in  the  area  depicted  in  Figure  3.3  was  treated 
as  a  random  variable  denoted  Y).  The  distribution  of  the  surface  heights 

represented  by  ^{X,  Y)  was  assumed  to  be  jointly  Gaussian.  This  is  a  valid  assump¬ 
tion  by  the  central  limit  theorem  as  long  as  the  value  of  n  is  kept  large  to  provide  a 
large  number  of  sample  points  (n^).  Ibr  the  central  limit  theorem  to  be  appropriate, 
the  sample  mean  and  variance  must  not  approach  extreme  values. 

Each  block  in  the  diaracterization  area  was  in  turn  treated  as  a  random  vari¬ 
able,  —  iSr,  Y  s  j6r)y  which  was  also  assumed  to  represent  a  joint  Gaussian 
distribution  of  heights.  K  the  Gaussian  assumption  was  not  made,  each  of  the  blocks 
in  Figure  3.3  would  have  to  be  sampled  at  increments  small  enough  to  allow  accurate 
development  of  the  governing  probability  distribution  £>r  each  block. 

This  procedure  would  become  almost  impossible  when  working  from  a  topo¬ 
graphic  map  because  of  the  various  sizes  of  the  homogeneous  areas  to  be  diaracter- 
ized.  The  size  of  the  homogeneous  areas  can  vary  from  several  square  IGlometers 
down  to  several  hundred  square  meters.  Unfortunatdiy,  the  contour  lines  are  only 
placed  at  three  meter  height  intervals  and  for  small  characterization  areas  the  distri¬ 
butions  for  the  individual  blocks  would  q>pear  as  uniform  distributions  at  a  single 
height.  EiXcept  fi>r  a  limited  number  of  cases,  this  would  probably  be  a  misleading 
characterization  of  the  height  distributions. 

One  can  extract  useful  information  about  the  properties  of  the 
tributions  based  on  the  assumption  that  they  are  all  Gaussian,  the  characteristics  of 
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Sr  -  Sr 


Ugoie  3.4.  Detenninuig  the  Sample  Point  Mean  hmm  Contonx  lines. 

the  information  provided  hy  the  contour  lines  the  topographic  map,  and  proper 
selection  of  the  spadng  Sr. 


S.S.l  Determining  the  Sample  Point  Mean  Value,  As  previously  mentioned, 
the  contour  lines  on  the  USGS  topographic  map  are  labeled  at  intervals  of  three  me- 
ters.  The  aze  of  the  areas  represented  by  the  Cihj)  di8tribnti|3ns  is  determined 
Vy  the  selection  of  the  quantity  Sr.  If  dr  is  selected  such  that  a  majority  of  these 
distribution  areas  coincide  ^thin  the  area  delimited  by  two  contour  lines,  a  rea¬ 
sonable  estimate  of  the  mean  value  of  the  individual  distributions 
figure  3.4  is  provided  to  illustrate  this  point. 

figure  3.4  illustrates  two  distribution  areas  labded  C(^i)  ' 

distribution  area  identified  as  C(*>i)  >>  bounded  by  two  oontour  fines ^ 


be  extracted. 

^:{i  +  lj).  The 
'labeled  Li  and 


Li+3m'  Rom  this  one  can  infer  that  a  majority  of  the  surface  heights  in  must 

lie  between  the  two  values  represented  by  the  contour  lines.  In  this  case,  the  mean 
value  of  C(t>j)  would  be  determined  by: 

1‘iJ  = (3-27) 

Even  if  great  care  is  taken  in  selecting  the  proper  value  of  6r,  there  will  still 
be  distribution  areas  which  engulf  more  than  two  contour  lines.  This  is  illustrated 
by  the  distribution  area  labeled  ^(i  + 1,  j)  in  Figure  3.4.  The  concept  described  by 
Equation  3.27  can  easily  be  extended  to  include  the  multiple  contour  line  case  as: 

+  i.  (3.28) 

Where: 

Lh  represents  the  contour  line  of  highest  value. 

Lt  represents  the  contour  line  of  lowest  value. 

Applying  Equation  3.28  to  the  distribution  area  in  Figure  3.4  which  is  labeled 
^(*  +  1, j)  would  yield: 

=  +  (3.29) 

By  carefully  selecting  the  proper  value  for  6r  and  using  the  procedure  outlined 
by  Equations  3.27  and  3.28  valid  mean  values  for  the  n*  blocks  in  each  site  can  be 
identified. 

S.S.2  Determining  the  Sample  Point  Variance.  The  variance  of  the  distribu¬ 
tion  areas  illustrated  in  Figure  3.4  can  also  be  determined  with  reasonable  accuracy 
from  the  data  provided  by  the  topographic  map  contour  Knes.  Using  the  distribution 
area  represented  by  C{i,j)  in  Figure  3.4,  one  can  again  infer  that  a  majority  of  the 
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saiface  heights  fall  between  the  values  identified  1^  the  oontonr  ]in«j  labeled  Li  and 
Li+3m-  For  a  Gaussian  distribution,  95.4  percent  of  the  distribution  values  fidl  within 
plus  or  minus  two  standard  deviations  of  the  mean  value.  To  derive  the  variance  of  a 
particular  area,  it  was  assumed  that  the  difference  in  contour  line  values  represented 
plus  or  minus  two  standard  deviations. 

4a  =  Li  —  Li+3n  (3.30) 

Solving  equation  3.30  £)r  a  and  noting  that  the  variance  equals  the  standard 
deviation  squared  jdelded  the  Allowing  formula  for  distribution  areas  bounded  by 
two  contour  fines: 


<rij  = 


_  (Ll+Sm  ~ 


16 


(3.31) 


Extending  Equation  3.31  to  the  cases  where  the  distribution  area  was  bounded 
fay  more  than  one  contour  fine  can  be  illustrated  ((t  +  IJ)  in  Figure  3.4.  fibr 
distribution  areas  similar  to  C(t  +  l«i))  the  plus  or  minus  two  standard  deviations  is 
now  represented  by  the  difference  in  height  of  the  highest  value  contour  fine  and  the 
lowest  value  contour  fine. 


^ij  = 


[Lk-L,) 

16 


(3.32) 


S.S.S  Estimating  the  Sample  Mean  and  Variance.  The  sample  mean  value 
was  estimated  as  the  average  value  of  the  mean  values  of  the  sample  point  dis* 
tributions.  The  estimate  of  the  sample  mean  was  determined  using  the  following 
equation: 


inn 


isl  jsl 


(3.33) 
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In  a  similar  manor  the  estimate  of  the  sample  vaii2mce  was  taken  as  the  average 
value  of  the  variances  of  the  individual  sample  point  distributions.  The  estimate  of 
the  sample  variance  was  determined  using  the  Allowing  equation: 


(3.34) 


S.S.4  Estimating  the  Sample  Correlation  Length.  The  oorrdation  length  is 
defined  as  the  distance  required  to  cause  the  correlation  function  reach  a  value  of 
The  correlation  function  fi>r  the  joint  Gaussian  distribution  as  it  iq>plies  to  the 
problem  at  hand  is  given  by: 


PaiJU(fn,n)  =  exp 


(3.35) 


ha  Equation  3.35  ({i,j),C{m,n)  are  any  two  of  the  n’ joint  Gaussian  distri¬ 
butions  which  are  part  of  C(-^,  Y).  The  quantity  /  represents  the  oorrdation  length 
or  distance  between  the  two  sdected  distributions.  The  quantity  r  represents  the 
distance  between  the  two  sdected  distributions  and  is  given  by: 


r^  =  {Xi-Xmy  +  {Yi-Y,)^ 


(3.36) 


The  correlation  function  for  any  two  distributions  can  also  be  expressed  as: 


PCijCmk 


^  [(C»i  ***  Pii)iCmk  *•*  Mmt)] 


(3.37) 


Where  in  Equation  3.37  the  £[*]  represents  the  average  or  expected  value  of 
the  argument.  Carrying  out  the  multiplication  of  terms  in  Equation  3.37  and  using 
the  &ct  that  the  expected  value  of  a  sum  is  equal  to  the  sum  of  the  expected  values, 
allows  the  correlation  function  to  be  expressed  as: 
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(3.38) 


For  any  given  combination  of  sample  distributions,  all  the  quantities  in  Equa¬ 
tion  3.38  can  be  determined  with  the  exception  of  E[CijCmk]’  Similarly  for  any  com¬ 
bination  of  of  surface  distributions,  the  only  unknown  quantity  in  Equation  3.37  is 
the  correlation  distance.  By  equating  the  two  equations  and  determining  an  estimate 
hi  ■^[Co'CmJb]  an  average  value  hi  the  correlation  distance  could  be  determined. 


During  the  literature  search,  the  researcher  could  find  no  accurate  method 
described  to  estimate  ^[CijCmfc]-  B  effectively  represents  Low  much  inter-dependence 
there  is  between  two  surface  heights  in  this  problem.  It  was  decided  to  estimate 
E[CiiU]  by  the  variance  estimator  for  the  sample  Making  this  substitution  into 
fk^uation  3.38  and  equating  Equations  3.37  and  3.38  yields: 


(3.39) 


Ihking  the  logarithm  of  both  rides  of  Equation  3.39  and  solving  for  produces 
the  following  expression: 


hi  f’lzHUhnk'] 

Equation  3.40  provides  a  value  for  the  cotrdation  distance  hr  a  single  combi¬ 
nation  of  surface  height  distributions.  However,  any  individual  distribution  could  be 
correlated  with  any  of  the  other  distributions.  Since  there  are  distributions,  there 
are  n*  possible  combinations  of  distributions.  'V^th  this  in  mind,  the  estimator  for 
the  average  correlation  distance  is  provided  below. 


3-21 


<1  n  n  n  n  ^ 

r^  =  ^EEEE-rr- — x 

”  i=i  j=i  m=i  *=i  ]n  ( S-znansihA 

Where  in  Equation  3.41  is  given  by  Equation  3.36.  A  FORTRAN  program 
titled  SiteChar.for  was  developed  to  determine  the  estimators  previously  described 
(/^,  and  /i)  horn  the  site  data  which  was  manually  extracted  from  the  topographic 
map. 

3.4  Determining  the  NRCS  of  the  Site  Characterization  Block. 

At  this  point  in  the  process,  the  topographic  map  has  been  annotated  at  each 
a£  the  selected  sites  with  the  rsmge  cell  Hnes  and  the  antenna  patterns  tor  both  the 
transmitter  and  receiver.  To  apply  the  concepts  outlined  in  the  previous  sections 
the  raw  data  must  be  extracted  from  the  map  at  each  site.  A  single  homogeneous 
area  at  each  site  was  selected  for  analysis.  IVom  each  cf  these  areas  three  range  cells 
were  selected. 

The  range  cells  selected  £br  site  one  were  cell  numbers  three,  £>ur,  and  five.  At 
site  two,  range  cell  numbers  five,  six,  and  seven  were  selected,  finally  at  site  three, 
range  cell  numbers  21,  22,  and  23  were  selected  ibr  the  calculations  that  follow. 

The  areas  to  be  characterized  were  too  small  to  allow  the  data  to  be  directly 
extracted  from  the  map.  Therefore,  the  selected  areas  of  the  three  map  sites  were 
enlarged  via  photocopier  to  allow  more  accurate  data  extraction.  The  three  areas 
received  different  degrees  of  enlargement  due  to  the  varying  density  of  the  contour 
fines  at  each  site. 

The  characterization  blocks  used  were  divided  into  400  equally  spaced  areas. 
The  value  of  400  was  selected  to  insure  enough  sample  points  were  used  to  insure 
the  assumption  that  the  random  variable,  ^{X,Y)y  was  jointly  Gaussian.  The  as¬ 
sumption  was  based  on  the  central  limit  theorem  which  requires  a  large  number  of 
samples.  A  more  realistic  picture  of  the  distribution  would  be  possible  with  even 
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inore  sample  points  however,  all  data  points  had  to  be  extracted  from  the  map  by 
hand.  The  value  of  400  sub  areas  was  a  compromise  to  minimize  the  manual  labor 
required  to  extract  data  at  each  site  and  still  maintain  the  validity  of  the  Gaussian 
assumption  which  was  based  on  the  central  limit  theorem. 

Site  one  was  divided  into  400  square  sub  areas  which  measured  23.348  m 
per  side.  The  total  area  encompassed  by  the  site  one  characterization  area  was 
218.052ii'm^.  Each  of  the  sub  areas  in  site  two  were  31.667  m  per  side.  The  total 
area  involved  for  site  two  was  401.12/f’m^.  Finally,  the  sub  areas  in  site  three  were 
37.703  m  per  side  for  a  total  area  of  SOS.Olifm^. 

The  input  data  to  the  FORTRAN  program  Sitechar.for  took  the  form  of  two 
real  numbers  per  sub  area.  The  first  number  represented  the  value  of  the  highest 
contour  line  bordering  the  area  and  the  second  number  represented  the  lowest  value 
contour  line  bordering  the  area.  With  these  inputs,  5jtec/iar./orcalculated  the  mean 
(p),  the  variance  (cr^),  and  the  correlation  length  (!)  estimates  for  the  sample.  Ta¬ 
ble  3.1  displays  the  results  of  running  the  program  Sitechar.for  on  the  three  sets  of 
input  data.  All  of  the  collected  and  calculated  data  is  provided  in  appendix  A. 

5.5  Calculating  the  Range  Cell  Power. 

The  clutter  power  from  each  range  cell  was  calculated  with  Equation  2.25  for 
each  range  cell  at  both  pol2Lrization3.  Before  Equation  2.25  could  be  utilized  the 
NRCS  for  each  site  characterization  area  had  to  be  determined.  A  second  FOR¬ 
TRAN  program  titled  Sigma.for  Was  written  to  determine  the  NRCS  from  the  site 
characterization  data  and  knowledge  of  the  locations  of  the  site,  transmitter,  and 
receiver. 

Once  the  NRCS  for  the  characterization  area  was  know,  the  clutter  block  areas 
were  determined  and  used  along  with  the  known  antenna  gains  to  determine  the 
clutter  power  from  each  block.  Summing  the  clutter  power  from  edl  blocks  located 
within  a  range  cell  yielded  the  total  clutter  power  from  the  range  cell  in  question. 
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Site  1 

Site  2 

Site  3 

Sr 

23.348m 

31.667m 

37.703m 

81.859m 

73.748m 

44.084m 

or* 

186.074m’ 

209.830m’ 

68.574m’ 

P 

45.687m’ 

82.365m’ 

116.587m’ 

Tkble  3.1.  Site  Characterization  Data. 
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Site  1 

Site  2 

Site  3 

2.539277 

3.262236 

0.0 

71.11379 

326.910 

0.0 

Tkble  3.2.  NRCS  Data  from  Sigma.for. 


The  output  from  the  program  Sigma.for  was  two  NRCS  values  for  the  charac¬ 
terization  area.  The  first  ht  both  transmitter  and  receiver  using  vertical  polarization 
and  the  second  ht  both  employing  horizontal  polarization.  The  input  files  used  and 
the  data  files  resulting  from  Sigma.for  were  included  in  Appendix  A.  Ikble  3.2  lists 
the  results  running  the  program. 

The  input  data  required  £>r  5tyma./or  included: 

1.  The  transmitter  to  clutter  cell  distance  {Dt). 

2.  The  receiver  to  dutter  odl  distance  (Dr). 

3.  The  sample  mean  estimate  (ft). 

4.  The  sample  variance  estimate  (a^). 

5.  The  relative  permittmty  of  the  dte  characterization  area  (er)*) 

6.  The  transmitter  ate  elevation  (h(). 

7.  The  receiver  site  devation  (hr). 

8.  The  correlation  length  of  the  sample  (/). 

All  three  of  the  characterization  sites  were  populated  by  trees  and  the  average 
value  for  the  relative  permittivity  was  taken  to  be  4.0  -I-  jO.45  (8:8).  The  site  areas 
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were  assumed  to  be  composed  of  non-magnetic  material  and  the  average  value  £)r 
the  relative  permeability  was  taken  to  be  1.0  +i0.0. 

S.5.1  Clutter  Power  for  Site  One.  The  dutter  power  for  three  range  cells  was 
determined  for  site  one.  The  selected  range  ndln  were  cells  three,  four,  and  five. 

When  using  vertical  polarization,  three  dutter  blocks  were  illuminated  by  the 
transmitter  and  receiver  antenna  patterns.  One  of  the  blocks  was  created  by  the 
intersection  of  the  two  antenna  main  beams.  The  remaining  two  blocks  were  created 
by  the  intersection  of  the  transmit  antenna  side  lobe  and  the  receive  antenna  main 
lobe.  Inserting  the  known  information  for  range  cell  three  into  Equation  2.25  yields: 

p  _  PA^GRm<r^ 

(4t)3 

Where: 

Cxm  =  The  transmit  antetma  main  lobe  gain  (39537.4). 

Gt*  =  The  transmit  antenna  side  lobe  gain  (34.43). 

G/im  =  The  receive  antenna  main  lobe  gain  (100.0). 

Aci  =  The  area  outlined  by  the  intersection  of  both  antenna  main  beams  and 
the  eqtiirange  contour  lines  which  define  range  odl  3  (200000m^). 

Ac2  =  The  first  area  outlined  by  the  transmit  antenna  side  lobe,  the  receive  an¬ 
tenna  main  lobe,  and  the  equirange  contour  lines  which  define  range  cell  3  (65390m^). 

Ac3  =  The  second  area  outlined  by  the  transmit  antenna  side  lobe,  the  receive 
antenna  main  lobe  and  equirange  contour  Hues  which  define  range  cell  3  (287500m^). 

Dti  Was  the  distance  from  the  center  of  the  first  dutter  Mock  (Aci)  to  the 
transmitter  (13450  m). 

Bt2  Was  the  distance  from  the  center  of  the  second  dutter  block  (Aci)  to  the 
transmitter  (12650  m). 


GrtnAel  GTMAe2  GtsAcS 


DIM 


D?M  Di,m 


a^r3 


(3.42) 
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Da  Was  the  distance  from  the  center  of  the  third  clatter  Uock  (i4c3)  to  the 
transmitter  (14425  m). 

Dri  Was  the  distance  from  the  center  of  the  first  dutter  block  (i4ci)  to  the 
receiver  (6550  m). 

Dr2  Was  the  distance  from  the  center  of  the  second  dutter  block  (>1^2)  to  the 
receiver  (7350  m). 

Drz  Was  the  distance  from  the  center  of  the  third  dutter  block  {Acz)  to  the 
receiver  (5575  m). 

<r"  =  The  NRCS  £ot  the  site  one  characterization  area  (71.11379). 

Substituting  these  values  into  Equation  3.42  and  performing  the  indicated  op¬ 
erations  yielded  a  received  dutter  power  of  30.209  mW  fi>r  range  cell  three  using 
vertical  polarization. 

ibr  the  case  of  both  receiver  and  receiver  using  horizontal  polarization,  range 
cell  three  contains  only  two  dutter  blocks.  The  first  was  defined  by  the  intersection 
of  the  antenna  main  beams  and  had  an  area  of  109375m’.  The  second  dutter  blodc 
was  defined  by  the  intersection  of  the  transmit  antenna  side  bbe  and  the  receive 
antenna  main  lobe  and  had  an  area  of  220000m’.  Using  these  value  in  Equation  3.42 
along  with  =  2.621544  yielded  a  dutter  power  of  5.900  mW. 

All  of  the  range  cells  fi>r  aU  three  sites  could  be  calculated  with  Equation  3.42. 
In  some  instances,  one  of  the  dutter  blocks  defined  by  the  intersection  of  the  transmit 
antenna  ade  bbe  and  the  receive  antenna  main  bbe  did  not  fall  within  a  particular 
range  cell.  The  tables  which  follow  summarize  the  pertinent  range,  area  measure¬ 
ments,  distances,  and  resulting  dutter  powers. 

Thble  3.3  provides  a  tabular  listing  of  the  dutter  block  areas,  transmitter 
to  target  and  receiver  to  target  distances  br  both  systems  operating  with  vertical 
polarization.  In  all  the  range  and  data  tables,  the  8]rmbols  should  be  interpreted  as 
defined  in  Equation  3.42. 
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Site  1  Vert.  Pol. 

Range  Cell  3 

Range  Cdl  4 

Range  Cell  5 

Aci 

200,000m* 

296,875m* 

68,750m* 

Dtt 

13,450m 

14,500m 

15,225m 

Drl 

6, 550m 

5,750m 

5,275m 

Aci 

65,390m* 

236,250m* 

150,938m* 

Dti 

12,650m 

13,450m 

14, 625m 

Pn 

7,350m 

6,800m 

5,875m 

Ac3 

287,500m* 

134,400m* 

0.0m* 

Da 

14,425m 

15,450m 

N/A 

Drl 

5, 575m 

4,800m 

- N/A 

Ikble  3.3.  Site  One  Aiea  and  Range  Data  hr  Vertical  Polarization. 


3-28 


Site  1  Vert.  Pol. 

Range  Cdl  3 

Range  Cdl  4 

Range  Cdl  5 

Ac\ 

29.210mlV 

48.409mlV 

12.082mlV 

Ac7 

7A6BfiW 

27.878/ilV 

20.181/xlV 

Ac3 

24.122/ilV 

0.0mlV 

Total 

29.261mTV 

48.461mW 

12.102mlV 

Table  3.4.  Site  One  Clutter  Block  Power,  Vertical  Polarization. 

Ibble  3.4  provides  the  calculated  dutter  Uock  poiwers  £>r  site  one  with  both 
transmitter  and  leceiver  operating  with  vertical  polarization. 

Ihble  3.5  provides  a  tabular  fisting  of  the  dutter  Uock  areas,  transmitter 
to  target  and  receiver  to  target  distances  ibr  both  systems  operating  with  vertical 
polarization,  hi  all  the  range  and  data  tables,  the  qrmbols  should  be  interpreted  as 
defined  in  Equation  3.42. 

Thble  3.6  provides  the  calculated  dutter  block  powers  hi  site  one  with  both 
transmitter  and  receiver  operating  with  horizontal  polarization. 

The  validity  of  these  range  cell  dutter  powers  is  unknown  because  the  sponsor 
was  unable  to  supply  actual  system  measurements  of  site  one. 
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Site  1  Horz.  Pol. 


Range  Cell  3 


Range  Cell  4 


Range  Cell  5 


Acx  109,375m2  2<)6,875m2  24,062m* 

Dll  13,450m  14,500m  15,225m 

Drt  6,550m  5,750m  5,275m 

Ac2  0.0m*  196,875m*  103,125m* 

Da  N/A  13,450m  14,625m 

Dr2  N/A  6,800m  5,875m 

Ac3  220,000m*  100,000m*  0.0m* 

Dt3  14,425m  15,450m  N/A 

Dr3  5,575m  4,800m  N/A 


Table  3.5.  Site  One  Area  and  Range  Data  for  Horizontal  Polarization. 


Site  1  Horz.  Pol. 

Range  Cdl  3 

Range  Cell  4 

Range  Cell  5 

Ac\ 

570.S87mW' 

1.729m  W' 

150.994/iW^ 

Ac2 

0.0miy 

829.537nW^ 

492.339niy 

Acz 

1.199^W 

640.867nW^ 

O.OmW'^ 

Tbtal 

571.586;«W^ 

1.730mW^ 

151.586/11^ 

Ikble  3.6.  Site  One  Qattei  Block  Power,  Horizontal  Polarization. 


S.5.2  Clutter  Power  for  Site  Two.  T^ble  3.7  provides  a  tabular  Hsting  of  the 
clutter  block  areas,  transmitter  to  target  and  receiver  to  target  distances  £3r  both 
systems  operating  vrith  vertical  polarization.  In  all  the  range  and  data  tables,  the 
symbols  should  be  interpreted  as  defined  in  Equation  3.42. 

Table  3.8  provides  the  calculated  clutter  block  powers  for  site  two  with  both 
transmitter  and  receiver  operating  with  vertical  polarization. 

Ihble  3.9  provides  a  tabular  fisting  of  the  clutter  block  areas,  transmitter 
to  target  and  receiver  to  target  distances  for  both  qrstems  operating  with  verticzil 
polarization.  In  all  the  range  and  data  tables,  the  symbols  should  be  interpreted  as 
defined  in  Equation  3.42. 

Table  3.10  provides  the  calculated  clutter  block  powers  fi>r  site  two  with  both 
transmitter  and  receiver  operating  with  horizontal  polarization. 

The  validity  of  these  range  cell  clutter  powers  is  unknown  because  the  sponsor 
was  unable  to  supply  actual  system  measurements  of  site  two. 

S.5.S  Clutter  Power  for  Site  Three.  Ihble  3.11  provides  a  tabular  fisting  of 
the  clutter  block  areas,  transmitter  to  target  and  receiver  to  target  distances  fbr 
both  systems  operating  with  vertical  polarization.  In  all  the  range  and  data  tables, 
the  symbols  should  be  interpreted  as  defined  in  Equation  3.42. 

Thble  3.12  provides  a  tabular  fisting  of  the  clutter  block  areas,  transmitter 
to  target  and  receiver  to  target  distances  fbr  both  ^sterns  operating  with  vertical 
polarization,  hi  all  the  range  and  data  tables,  the  symbols  should  be  interpreted  as 
defined  in  Equation  3.42. 

Thbles  3.11  and  3.12  were  included  for  completeness.  The  power  level  calcula¬ 
tions  fsr  site  three  could  not  be  performed  because  the  expression  fbr  the  NRCS  at 
sight  three  could  not  be  determined. 
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Site  2  Vert.  Pol. 

Range  Cdl  5 

Range  Cdl  6 

Range  Cell  7 

Aci 

210,000m* 

245,000m* 

186,046m* 

Dn 

9, 87Sm 

10, 925m 

11,875m 

Drl 

10,625m 

9, 825m 

8, 875m 

Aoi 

121,875m* 

132,187.5m* 

131,250m* 

Dt2 

9,425m 

10,225m 

11,175m 

Dr2 

11,075m 

10,525m 

9, 825m 

Ac3 

132,812.5m* 

134,922m* 

65,609m* 

Da 

10,675m 

11,725m 

12, 600m 

1 

1 

9,825m 

9, 025m 

8,400m 

Tkble  3.7.  Site  Two  Area  and  Range  Data  fcr  Vertical  Polarization. 
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site  2  Vert.  Pol. 


Range  Cell  5 


Range  Cell  6 


Range  Cell  7 


Act 

99.400miy 

110.804mW 

87.280m 

Ac7 

Hi 

■ 

51.790/xW 

49.404AiVr 

Acs 

54.785^1^ 

54.675/iW 

26.576/iW^ 

Total 

110.911mW  I 

1 

87.356m 

Tkble  3.8.  Site  Two  CHutter  Block  Power,  Vertica 

Polarization. 
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site  2  Horz.  Pol. 

Range  Cell  5 

Range  Cell  6 

Range  Cell  7 

Aci 

159,687.5m* 

245,000m* 

163,750m* 

Dn 

9,875m 

10,925m 

11,875m 

Drl 

10,625m 

9, 825m 

8, 875m 

Ac2 

46,094m* 

132, 187.5m* 

131,250m* 

Di2 

9,425m 

10,225m 

11, 175m 

Dt2 

11,075m 

10,525m 

9, 825m 

Ac3 

132,812.5m* 

134,922m* 

37,500m.* 

Dtz 

10,675m 

11,725m 

12, 600m 

DrZ 

9,825r" 

9, 025m 

8, 400m 

Table  3.9.  Site  Two  Area  and  Range  Data  for  Horizontal  Polarization. 
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site  2  Horz.  Pol. 

Range  Cell  5 

Range  Cdl  6 

Range  Cell  7 

Aci 

754.263/ilV 

l.lOOmW’ 

799.994AiH^ 

Ac7 

191.560nW^ 

516.812nir 

'193.006nW^ 

Acs 

546.702nlV 

545.601nW^ 

151.581nW^ 

Total 

755.001/iir 

1.107mW^ 

T^ble  3.10.  Site  Two  Clutter  Block  Power,  Horizontal  Polarization. 


Site  3  Vert.  Pol. 

Range  Cell  21 

Range  Cell  22 

Range  Cell  23 

Ac\ 

97,500m’ 

97,500m’ 

97,500m’ 

Dt, 

15,325m 

15,625m 

15,950m 

Drl 

9, 925m 

9, 850m 

9,800m 

Ac2 

35,000m’ 

64,471.4m’ 

65,000m’ 

Dti 

14,825m 

15,150m 

15,450m 

Drt 

10,425m 

10,350m 

10,300m 

Acz 

69,062.5m’ 

69,062.5m’ 

69,062.5m’ 

Da 

15, 850m 

16, 150m 

16,475m 

DrZ 

9,400m 

9,350m 

9,275m 

Ikble  3.11.  Site  Three  Area  and  Range  Data  ibr  Vertical  Fol^Lrization. 
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Site  3  Horz.  Pol.  Range  Cell  21 


Range  Cell  23 


97,500m^ 


97,500m’ 


97,500m’ 


Dti  15, 325m  15, 625m  15, 950m 

Dri  9,925m  9,850m  9,800m 

Ac2  937.5m’  52,500m’  65,000m’ 

Dt2  14,825m  15,150m  15,450m 

Dri  10,425m  10,350m  10,300m 

Ac3  69,062.5m’  69,062.5m’  69,062.5m’ 

Dt3  15,850m  16,150m  16,475m 

Drz  9,400m  9,350m  9,275m 


Ikble  3.12.  Site  Three  Area  and  Range  Data  for  Horizontal  Polarization. 


IV.  Conclusions  and  Recommendations 


.f.l  Conclxissiona 

It  was  impossible  to  establish  how  accurate  the  process  proposed  by  this  re¬ 
search  was  due  to  the  lack  of  actual  data  as  a  source  of  comparison.  However,  the 
power  levels  calculated  for  the  selected  range  cells  at  sites  one  and  two  were  not  en¬ 
tirely  tmreasonable  for  the  v-utter  power  and  distances  involved.  This  in  itself 
was  a  degree  of  success. 

The  research  did  providf;  some  insight  into  the  complexity  of  realistically  im¬ 
plementing  the  process  and  iodicated  areas  requiring  fiirther  research.  The  most 
obvious  of  these  areas  was  discovered  with  the  site  three  data.  Barrick’s  very  rough 
surface  formulation  fails  for  the  special  case  when  dealing  with  the  combination  of 
low  grazing  angles  and  a  bistatic  angle  which  approaches  90”. 

Recall  from  lOiapter  2  the  NRCS  was  determined  from  the  expression: 

=1  M'J  (4.1) 

Where  in  Equation  4.1,  represents  a  modified  scattering  matrix  dement 
and  the  J  term  the  fidd  component.  For  the  Joint  Gaui^sian  distribution,  J  was 
given  by: 


(4.2) 


Where  the  quantities  in  Equation  4.2  s’,  (y,  and  have  been  defined  as: 


^  =  sin  0,-  —  sin  cos 


(4.3) 
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(y  =  sin  6i  sin  (ft. 


(4.4) 


—  COS  Oi  —  008  (4.5) 

At  low  grazing  angles,  (  dy  <  5°  where  Og  has  been  defined  as  f  —  0i  ),  the 
exponential  argument  blossoms  rapidly.  This  is  due  to  the  in  the  denominator 
of  the  exponent’s  argument.  For  small  values  of  6i  Og  approaches  |  and  hence  the 
cosine  oi  Og  t  .  oaches  zero.  Similarly,  as  the  bistatic  angle  approadies  |  so  does 
Of,  When  th**  e  two  situations  occur  simultaneously,  rapidly  approaches  zero  and 
renders  the  calculation  of  questionable  value. 

A  second  complication  with  the  process  was  the  actual  data  extraction  from 
the  topographic  map.  The  process  was  extremely  time  consuming  when  performed 
manually,  fbr  the  three  sites  in  the  project,  it  took  i^proximately  twelve  hours  to 
extract  the  data  and  enter  it  into  a  data  base.  Fbr  the  process  to  realistically  be 
used,  the  entire  target  area  would  have  to  be  characterized  which  in  turn  implies 
extracting  data  from  tho  map  fr>r  the  entire  radar  coverage  area. 

The  process  vras  unable  to  include  the  effects  of  small  scale  toughness  on  the 
received  clutter  power  due  to  the  lack  of  detailed  information  on  the  terrain  popula¬ 
tion.  It  is  unknown  what  the  magnitude  of  this  contribution  could  amount  to.  If  it 
is  not  a  significant  contributor  to  the  over  all  range  cell  dutter  power,  the  costs  of 
collecting  adequate  data  to  establish  its  level  may  fitr  exceed  the  return. 

Finally,  at  low  grazing  angles,  a  shadowing  function  is  almost  imperative.  The 
problem  again  is  that  the  shadowing  function  is  of  little  use  without  knowledge  of 
all  the  terrain  between  the  transmitter,  receiver,  and  target. 
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4.2  Recommtndations 

The  first  priority  would  be  to  establish  the  relative  validity  of  the  dutter  powers 
by  direct  oomparison  with  measured  values.  This  would  be  a  much  more  efRdent 
process  if  the  hardware  were  located  AFIT  where  the  researched  has  access  to  it. 

Manually  extracting  topographic  data  has  a  prohibitive  man  hour  cost.  The 
collection  of  height  data  for  the  very  rough  surface  scattering  must  be  automated. 
One  possible  approach  to  this  end  would  be  the  use  of  DTED  data.  However,  the 
digitized  data  comes  in  a  half  inch  tape  format  and  requires  a  relatively  long  lead 
time  to  acquire. 

If  digital  site  data  were  available,  the  homogeneous  terrain  areas  could  be  used 
to  establish  thdr  probability  distribution  on  an  area  by  area  basis.  This  would 
also  require  new  solutions  to  the  radiation  integral  in  Chapter  2  which  establishes  a 
value  ht  J.  The  automated  data  could  also  be  useful  £>r  establishing  the  validity 
of  the  Gaussian  height  distribution  assumption  which  was  based  solely  on  the  loose 
interpretation  of  the  central  limit  theorem. 
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Appendix  A.  Software 


The  following  FORTRAN  program  titled  Sitechar.for,  was  used  to  calculate 
the  sample  mean,  variance,  and  correlation  length  estimators. 

The  input  data  files  consisted  of  sets  of  numbers  vdiich  represented  the 
highest  contour  line  value  and  the  lowest  contour  line  value  which  delimited  each 
sub  area  in  the  sample.  The  three  input  data  files  ( Sitel.dai,  Site2.dat,  and  Site3.dat 
)  are  included  in  Appendix  B.  The  results  of  the  computations  are  included  in  three 
data  files  (  Sitelprb.dat,  Sitc2prb.dat,  and  Site3prb.dat }  which  are  also  qrovided  in 
Appendix  B. 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

c  SITECHAR.FOR 

C 

c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


IMPLICIT  NONE 

variable  declarations 

integer 

i 

!  loop  index 

integer 

j 

!  loop  index 

integer 

n 

!  loop  index 

integer 

m 

!  loop  index 

teal*8 

logarg 

realms 

high 

!  upper  contour  line 

A-1 


teal>»8 

low 

!  lower  contour  line 

ieal4>8 

nsamp 

!  number  of  samples  per  axis 

ieal*8 

mn 

I  mean  estimator 

teali»8 

isqrd 

1  distance  between  dist 

real*8 

hsqrd 

!  variance  estimator 

real«8 

kqrd 

!  correlation  length  est 

ieal*8 

ddr 

I  sample  separation  distanc 

ieal>»8 

aeta(500,500) 

!  mn  Z(i  j) 

ieal>t<8 

agma(500,500) !  standard  deviation  2^iJ) 

ieal<i>8 

!  temp  storage  xeg 

character*  15 

wifilnmn 

!  file  with  height  samps 

character*  15 

outfilnam 

!  file  with  results 

initialize  uatiables 

iisamp=:20.0 

delr=:10.0  ' 

li8qrd=:0.0 

xnii=0.0 

lfqid=0.0 

get  input  data  file  name 

'write(*,440}*  Enter  input  data  file 
read(*,450)  infilnam 


I  Kt  to  number  of  samples 
!  sample  sep  distance 


get  output  file  name 


c 

c 

c 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


Write(+,440)’  Enter  output  data  file  name:’ 
read(*,450)  outfilnam 

get  number  divisions  per  axis 

write(*,440)’  Enter  number  of  divisions  per  axis:’ 
read(*,480)nsam.p 

get  sample  w’dth 


write(*,440)’  Enter  sample  width  :’ 
read(*,400)  delr 


open  input  data  file 


open(unit=38,file=infilnam,status=’OLD’) 


open  output  file 

open(unit=39,file=outfilnam,status=’NEW’) 


read  in  the  contour  line  values  and  echo  to  screen 


do  100  i=l,nsamp,l 
do  110  j=l,nsamp,l 
read(38,err=1000,£mt=480)  high 
read(38,err=1000,£mt=480)  low 
zeta(ij)=((high-Iow)/2.0)+low 


!  set  up  the  row  counter 
!  set  up  the  col  counter 
!  highest  contour  line 
!  lowest  contour  fine 
!  find  mu(ij)  and  save 
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agma(ij)=(lugh-low)/4.0  !  find  8igina(iJ)  and  save 

C  wxite(4>,400)  zeta(ij) 

C  write(*,400)  sigma(ij) 

110  continue 
100  continue 
C 

C  calculate  the  mean  estimate 

C 

mu=0.0 

do  120  i=l,nsamp,l 
do  130  jsl^amp,! 
mu=mu+zeta(ij) 

130  continue 
120  continue 
mu=mu/(nsamp<i»t'2) 
write(39,410)  ’mu  =  ’,mu 
wiite(*,410)  ’mu  =  ’,mu 
C 

C  calculate  the  height  variance  estimate 

C 

hsqrd=0.0 

do  140  isl,nsamp,l 

do  150  j=:l^amp,l 

h8qtd=hsqrd+(zeta(ij)-mn)>t»t>2 

150  continue 

140  continue 

hsqrd=h8qrd/(nsamp**2) 
write(39,420)  ’hsqrd  ’^isqrd 
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wtite(4<,420)  lisqrd  =  ’,lisqrd 
calculate  the  correlation  length  estimate 
lsqrd=:0.0 

do  160  i=l,n8amp,l 
do  170  j=l^amp,l 
treg=0.0 

do  180  m=l,n8amp,l 

do  190  n=l^amp,l 

r8qrd=((i-m)*delr)**24-((j-*i)*<lelr)**2 

logarg=(hsqrd-zeta(ij)%zeta(m,n)) 

logarg=ABS(logarg/(sigma(iJ)>Ksigma(m^))) 

treg=r8ord/LOG(Iogarg) 

190  continue 
180  continue 

l8qrd=:lsqrd+treg/(n8amp*4>2) 

170  continue 
160  continue 
lsqrd=Isqrd/ (nsamp4>«2) 

«xite(39,420)  Isqrd  =  ’^sqrd 

Save  the  sample  means 

do  200  i=l,nsamp,l 
do  210  j=l^amp,l 
write(39,490)  ’mu(’4,’,’j,’)  =  ’,zeta(ij) 

210  continue 


200  continue 


Save  the  sample  standard  deviations 


do  220  i=l,nsamp,l 
do  230  j=l^amp,l 

vnite(39,500)  ’sigma(’4,’,’ j,’)  =  ’,sigma(i  j) 
230  continue 
220  continue 
goto  999 


list  of  faimat  statements 

400  fi>rmat(F8.3) 

410  fi>rmat(A6,F8.3) 

420  £>nnat(A10,F8.3) 

430  £)rmat(A25) 

440  fi}rmat(A35) 

450  £)rmat(A15) 

460  £>rmat(I4) 

470  fi)nnat(A10,F15.3) 

480  fi>nnat(F3.0) 

490  fi>rmat(A342,A142,A4,F8.3) 
500  £)rmat(A642,A142,A4,F8.3) 


detected  ecroi  actions 


1000  imte(*,430)M!!  DATA  READ  ERROR  !!!’ 


c 

c 

c 


({oto  999 

1010  write(*,430)M!!  DATA  V/RITE  ERROR  !!!’ 


dose  the  data  files  and  exit  the  program 

999  dose(unit=38) 
do8e(unit=39) 
write(*,410)  ’mu  =  ’,mu 
write(*,420)  ’hsqrd  =  ’,hsqrd 
wiite(*,420)  ’Isqrd  =  ’4sqrd 
stop 
end 

The  program  implements  Barrick’s  rough  surface  scattering  fermulas  &om  the 
stistics  of  the  characterization  area  generated  with  5ttechar./orand  the  transmitter, 
receiver,  and  target  area  geometry. 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

c 

C  SIGMA.FOR 

C 
C 
C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

IMPLICIT  NONE 

c 

C  variable  deduations 
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/ 


c 


integer 

i 

!  bop  index 

integer 

n 

!  loop  index 

integer 

m 

!  loop  index 

complex 

ctem 

real*8 

j 

! 

real+8 

al 

! 

real=i'8 

a2 

! 

reu!+8 

a3 

! 

rw'il+8 

a4 

! 

real*8 

Isqrd 

!  correlation  kngth 

real*8 

angi 

! 

real*8 

treg 

!  temp  storage  reg 

real+8 

thetai 

!  incident  wave  ang 

real*8 

thetas 

!  scattered  wave  an 

reai*8 

phis 

!  scattered  wave  ph 

real*  8 

Dt 

!  linear  dist  trans-c! 

real*8 

Dr 

!  linear  dist  rcvr-cl 

reali'8 

Db 

!  linear  disi  trans-r 

real*  8 

ht 

!  transmitter  height 

real*  8 

hr 

!  receiver  height 

real*8 

he 

!  clutter  mean  plan 

real*8 

sigmaw 

real*8 

si^nahh 

1 

complex 

mux 

!  relative  permeabil 

complex 

epsilonr 

!  relative  permittivi 

complex 

Rper 

complex 

Rpar 
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ieal*8 

betaw 

ieal*8 

betahh 

real*8 

ssqrd 

! 

teal*8 

zetax 

1 

real«8 

zetay 

f 

teal*8 

I» 

teal>t<8 

cti 

!  COS(thetai) 

realt8 

sti 

!  SIN(thetai) 

real*8 

cts 

!  COS(thetas) 

ieel*8 

sts 

!  SIN(thetas) 

Teal>t>8 

cps 

!  COS(phis) 

Tealv8 

sps 

!  SIN(phis) 

complex 

rooti 

zeal*8 

zetaz 

real4<8 

si 

!  SIN(angi) 

rsal*8 

d 

!  COS(angi) 

zeal*8 

nsamp 

!  number  of  dutter 

ieal*8 

beta 

teal*8 

deltar 

?  width  of  dut  cell 

zeal*8 

hsqrd 

!  variance  of  dutter 

character  15 

iTifilnii.Tn 

!  file  with  height  sa 

character*  15 

outfilnam 

!  file  with  results 

initialize  variables 

Iit=97.572 

lu-144.0 

Db=19248.0 


!  set  transmitter  he 
!  set  receiver  height 
!  set  trans-rcvr  dist 


pi=3.141592653589793 


get  input  data  file  name 

iRite(*,440)*  Enter  input  data  file  name  :* 
read(*,450)  infitnam 

get  output  file  name 

wiite(«,440)*  Enter  output  data  file  name:’ 
iead(*,450)  outfilnzun 

open  input  data  file 

open(unit=:38^e=infilnam^tatttsss’OLD’) 
open  output  file 

open(unit=39^e=outfilnam^tatus=:  ’NEW’) 

read  in  transmitter  elevation  and  echo  to  screen 

xead(38,err=1000,fi3it=400)  Dt 
write(39.480)’Dt  =  ’,Dt 
write(*,480)’Dt  =  ’,Dt 

read  in  receiver  devation  and  echo  to  screen 


zead(38,err=1000,fmt=400)  Dr 
wtite(39,480)’Dr  =  ’,Dt 
VTrite(*,480)’Dr  =  ’,Dr 

lead  in  mean  duttar  height  aud  echo  to  screen 

read(38.,e:r=1000,{int=400)  he 
write(39,480)’hc  =  ’^c 
write(*,480)Tic  =  ’4»c 

read  in  relative  permeability  and  echo  to  screen 

2ead(38,err=1000,bnt=490)  mur 
write(39,500)’mn  =  ’,mux 
write(*,5(K))’mn  =  ’^ur 

read,  in  relative  permittivity  and  echo  to  screen 

read(38,err=1000,fmt=490)  epsilonr 
write(39,500)*ep8  =  epsilonr 
write(*,500)’eps  —  epsilonr 

read  in  area  of  characterization  site  and  echo  to  screen 

tead(38,err=1009,£mt=4C0)  nsamp 
Viite(39,480)’DIV  =  nsamp 
write(*,480)’EIV  =  nsamp 
read(38,err=1000,fmt=400)  deltar 
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write(39,480)’DEL  =  *,ddiar 
write(*,480)’D£L  =  ',<ldtar 

read  in  clutter  variance  and  echo  to  screen 

zead(38,err=1000,£ait=400)  hsqrd 
'write(39,480)’h8qrd  =  *, hsqrd 
'write(*,480)’h8qrd  =  *4>*qrd 

read  in  correlation  distance  and  echo  to  screen 

iead(38,err=1000,&it=400)  Isqrd 
irrite(39,470)’lsqrd  = 
write(*,470)lsqrd  = 

calculate  the  scattering  angles 

thetai=(pi/2.0)-DATAN((ht.hc)/Dt) 
thetas=(pi/2.0)-DATAN((hr-hc)/Dr) 
beta=DACOS((Dt**2+Di|>»*2-Db*v2)/(2.0*Dt*Dr)) 
phis=pi-beta  1 

vnite(39,470)  ’tL  'tai  =  ’,thetai 
write(39,470)  ’thetas  =  thetas 
vrrite(39,470)  ’beta  =  ’,beta\ 
vrrite(39,470)  ’phis  =:  ’,phis  1 
write(«,470)  ’thetai  =:  ’,thetali 
write(*,470)  thetas  =  *, thetas 
write(*,470)  ’beta  =  *,beta 
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c 

c 

c 


c 

c 

c 


c 

c 

c 


write(*,470)  ^his  =  ’tphis 
determine  the  angle  i 

tteg=1.0-DSIN(thetai)*DSIN(thetas)*DCOS(pIus) 

treg=(treg+DCOS(thetai)*DCOS(theta8))/2.0 

angi=DACOS(DSQRT(treg)) 

write(*,470)  ’angle  i  =  *,angi 

write(3d,470)  ’angle  i  =  ’,angi 

determine  the  a  ooefEdents 

cti— DCOS(thetai) 

8U=DSIN(thetai) 

cts:i:DCOS(thetas) 

sta=DSIN(theta«) 

q>8=:DCOS(phi8) 

^8=DSIN(phi8) 

al=1.0+8ti*8t8*q)8-cti4>ct8 

a2=Cti*8t8+8ti4<Ct84>q)8 

a3=8ti<»ct8+cti*8t8*q)8 

a4=cti+ct8 

I 

determine  the  reflection  ooefficient8 


d=DCOS(angi) 

8i=:DSIN(angi) 

ctemp=:(1.0,0.0)>*8i 


A43 


iooti=CSQRT(mur*ep8ilonr-(ctemp*4>2)) 

iper=(mttr*ci-rooti)/(miir4'ci+rooti) 

ipar=(epsiloiii4>ci-rooii)/(ep8iloni*ci+tooti) 

determine  aeta  x,7,  and  % 

Ktaz=8ti-8t8>t>q>s 

aeta7=st8*8p8 

aeias=0.0-cti*q>s 

determine  s*<t>2 

■qrd=4.0*li8qrd/l8qrd 
write(39,470)’88qrd  s  ’,88qrd 
wxite(*,470)’88qrd  =  *,8sqrd 

determine  J 

ja:(DI2XP(0.0-(zetax’i>*2+aeta7*<i'2)/(ssqrd>t>zetaz*4>2))) 

j=:j*4.0/(Mqrd>»zetaz>«>4>2) 

wiitc(39,480)’j  =  ’j 

wtite(*,480)’j  = 

determine  the  scattering  matrix  components 

ctemp=(a2*a34>rpar+sti*sts*(sps**2)*rper)/(al*a4) 

betavT=CABS(ctemp) 

Ctemp=(0.0<sti4>st8*(sps4>4>2)*rpar>a2*a3*rper)/(al*a4) 


betahli=CABS(ctemp) 
write(39,470)’b€taw  =  ’,betaw 
imte(*,470)’betaTT  =  ’,betaw 
Trrite(39,470)’betabb  =:  ’,betabb 
write(4>,470)’betahh  =  ’,betahb 


determine  the  NRCS 

agmaTV=(betaTY>i"«>2)4>j 
■igmahh=  (betidih<»*2)  *j 
write(39,470)’8igmaTy  =  *,sigmavT 
imte(*,470)*rigmaTT  =  *,ngmaTT 
mite(39,470)’sigmahh  =:  ’^'gmahh 
t7rite(*,470)*8igmalLh  =  *^gmahh 

list  of  format  statements 


400  format(F10.3) 

410  format(A6,Fl|).3) 

420  format(A10,F10.3) 
430  format(A25) 

440  format(A30) 

450  format(A15) 

460  format(I3) 

470  formate  ’,A10,Ejl3.7) 
480  formate  ’,A9,E13.7) 
490  fonnat(F10.5^10.5) 


500  fonnat(A9,Pl2.7,F12.7) 

C 

C  detected  error  actions 

C 

1000  wiite(*,43U)M!!  DATA  READ  ERROR  !!!' 
1010  wzite(*,430)M!!  DATA  WRITE  ERROR  !!!’ 
C 

C  dose  the  data  files  and  exit  the  program 

C 

999  dose(anit=38) 

dose(anit=39) 

stop 

end 
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Appendix  B.  Topographic  Data 

The  data  which  Ssllows  u  the  ate  one  contour  line  information  used  to  estimate 
the  ate  characterization  area  statistics. 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCrCCC' '( 


7=1 

7=2 

7=3 

y=4 

7=5 

Z=1 

60  57 

60  57 

60  57 

63  57 

63  57 

x=2 

60  57 

60  57 

63  57 

63  57 

63  57 

x=3 

63  60 

63  60 

69  60 

69  60 

69  60 

66  60 

66  60 

66  63 

66  63 

72  63 

x=5 

66  60 

66  63 

69  63 

69  63 

69  66 

x=6 

66  63 

66  63 

69  63 

6S  66 

72  66 

x=7 

66  63 

69  63 

69  66 

72  66 

75  69 

x*8 

69  60 

72  66 

72  66 

75  69 

75  69 

X=:9 

72  66 

75  69 

78  69 

78  72 

81  72 

X=10 

75  69 

78  72 

81  75 

84  75 

84  78 

X=:ll 

75  69 

81  75 

84  78 

84  81 

87  81 

X=12 

75  72 

81  75 

84  78 

87  81 

90  84 

X=13 

75  72 

81  75 

84  78 

90  81 

93  84 

x=14 

75  72 

81  75 

84  81 

90  81 

96  87 

x=15 

78  72 

84  75 

87  81 

93  84 

96  90 

x=13 

78  72 

84  78 

90  81 

93  87 

96  90 

x=17 

81  75 

84  81 

90  84 

93  87 

96  90 

x=18 

81  75 

84  78 

87  81 

90  84 

93  87 

z=19 

81  75 

84  78 

84  81 

37  81 

87  81 

x=20 

78  72 

81  75 

81  75 

81  75 

81  75 
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cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


y=6 

y=7 

y-3 

y=9 

y=10 

X=1 

63  CO 

63  60 

66  CO 

66  63 

66  63 

x=2 

63  60 

66  69 

66  63 

69  63 

69  66 

x=3 

69  CO 

72  63 

72  63 

72  66 

72  66 

x=4 

72  63 

69  66 

69  uC) 

72  63 

72  69 

x=5 

72  66 

72  66 

72  66 

75  66 

75  69 

x=6 

72  69 

75  69 

75  69 

75  72 

78  72 

x=7 

75  69 

75  69 

78  69 

78  72 

81  75 

X 

II 

CO 

75  69 

78  72 

78  72 

78  75 

81  75 

x=9 

81  75 

81  75 

84  75 

87  78 

90  78 

x=10 

84  78 

87  78 

90  81 

93  84 

102  87 

II 

X 

90  81 

93  84 

96  87 

102  90 

105  93 

II 

K 

93  84 

99  90 

102  90 

108  99 

111  102 

x=13 

99  90 

105  93 

108  99 

111  105 

114  108 

11 

K 

102  90 

105  99 

111  102 

114  108 

114  111 

II 

K 

105  96 

108  102 

111  105 

111  105 

111  105 

to 

II 

H 

102  96 

105  99 

111  102 

105  99 

105  99 

x=17 

102  93 

102  93 

102  90 

102  30 

99  87 

00 

II 

K 

96  87 

96  87 

93  84 

90  81 

90  78 

x=19 

87  81 

87  78 

84  78 

84  75 

81  75 

O 

11 

H 

81  75 

81  72 

78  72 

78  75 

81  75 
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/ 


y=ll 

y=12 

y=13 

7=14 

7=15 

Z=1 

66  63 

66  63 

66  63 

66  63 

66  63 

x=2 

69  63 

69  63 

69  63 

72  63 

72  63 

x=3 

72  66 

69  66 

72  66 

72  69 

72  69 

x=4 

75  69 

75  69 

78  69 

78  72 

78  72 

3Ps5 

75  72 

78  72 

78  72 

78  75 

78  75 

x=6 

81  72 

81  75 

81  75 

84  75 

84  75 

x=7 

84  75 

87  78 

87  81 

90  81 

90  84 

x=8 

84  75 

90  78 

93  81 

96  84 

96  87 

x=9 

96  84 

102  90 

102  90 

102  93 

102  96 

x=10 

105  90 

105  99 

108  99 

108  99 

108  102 

z=ll 

108  102 

111  105 

111  105 

111  105 

108  105 

x=12 

114  105 

114  108 

114  108 

111  108 

111  105 

xsl3 

114  111 

114  111 

114  108 

111  105 

108  102 

X=:14 

114  111 

114  108 

114  111 

114  108 

114  108 

x=15 

111  102 

111  105 

111  102 

108  96 

102  87 

z=16 

105  96 

105  93 

102  90 

96  87 

93  84 

x=17 

96  84 

96  84 

93  84 

90  84 

87  81 

x=18 

87  81 

87  81 

87  81 

84  78 

84  81 

x=19 

84  81 

84  81 

84  81 

84  78 

84  78 

x=20 

84  78 

81  78 

84  75 

81  75 

81  75 
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r=i6 

p=17 

y=18 

7=19 

7=20 

Z=1 

66  63 

66  63 

66  63 

66  63 

66  63 

x=2 

69  63 

69  63 

66  63 

69  63 

69  63 

x=:3 

72  66 

72  66 

69  63 

69  66 

69  66 

x=4 

78  69 

75  69 

75  66 

72  66 

72  66 

x=5 

78  75 

81  72 

81  69 

75  69 

75  69 

x=6 

84  75 

84  78 

84  78 

81  75 

81  72 

x=7 

87  81 

87  81 

90  81 

87  78 

84  73 

x=8 

99  87 

99  87 

96  87 

93  84 

87  81 

x=9 

102  96 

99  96 

99  90 

96  87 

90  84 

x=10 

105  99 

102  96 

99  93 

96  87 

90  81 

x=ll 

105  102 

102  99 

99  93 

93  87 

87  81 

x=12 

105  99 

102  96 

99  90 

93  84 

87  81 

Xal3 

105  99 

99  93 

96  87 

90  84 

87  81 

x=14 

108  102 

102  90 

93  81 

87  81 

84  78 

x=15 

96  84 

90  81 

87  81 

84  78 

81  72 

90  81 

84  81 

84  78 

81  78 

81  75 

x=17 

84  81 

84  81 

84  78 

81  75 

78  72 

xsl8 

84  78 

84  78 

81  75 

78  72 

75  72 

Z=:19 

81  75 

78  72 

78  72 

75  72 

75  69 

x=20 

81  75 

75  72 

75  69 

75  69 

75  66 
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Tlie  data  which  follows  is  the  site  two  contour  line  informatioa  used  to  estimate 
the  die  characterization  area  statistics. 
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y=l 

y=2 

II 

y=4 

li 

X=1 

87  C9 

90  78 

93  81 

93  84 

93  84 

II 

to 

93  78 

96  87 

96  90 

96  90 

96  90 

x=3 

93  78 

99  93 

102  93 

102  96 

99  93 

w 

II 

H 

99  87 

102  96 

102  99 

102  96 

99  96 

x=5 

105  96 

105  99 

105  99 

102  99 

99  96 

x=6 

105  99 

105  102 

105  99 

102  99 

99  96 

x=7 

108  102 

103  102 

105  102 

105  99 

102  96 

x=8 

108  102 

108  102 

105  102 

105  99 

102  96 

o 

II 

K 

108  99 

108  102 

105  99 

102  96 

102  93 

o 

II 

H 

102  90 

102  96 

102  96 

102  93 

99  90 

x=ll 

96  87 

99  93 

99  93 

99  90 

96  87 

X 

II 

to 

93  81 

96  90 

96  90 

96  87 

90  81 

x=13 

90  81 

93  87 

93  84 

90  81 

87  78 

x=14 

90  78 

90  78 

90  75 

84  72 

81  72 

x=15 

78  72 

78  72 

78  59 

75  69 

75  69 

x=16 

72  63 

72  63 

72  C6 

72  66 

72  63 

x=17 

69  GO 

66  CO 

66  CO 

66  60 

66  60 

x=ri 

63  CO 

63  GO 

63  57 

63  57 

63  57 

x=19 

60  54 

CO  5  4 

CO  57 

63  CO 

60  57 

x=20 

CO  57 

CO  57 

63  57 

63  60 

63  CO 
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ys=6 

7=7 

y=8 

y=9 

7=10 

Z=1 

90  81 

87  78 

84  75 

84  75 

81  72 

x=2 

93  87 

93  84 

90  81 

87  78 

84  78 

x=3 

96  90 

96  90 

93  87 

90  84 

90  81 

x=4 

99  93 

96  93 

96  90 

93  87 

90  87 

x=5 

99  93 

96  93 

96  90 

93  90 

90  87 

r=6 

99  93 

96  93 

96  93 

93  90 

93  87 

x=7 

99  96 

99  93 

96  93 

93  90 

93  87 

x=8 

99  96 

99  93 

96  90 

93  87 

90  84 

x=9 

99  90 

96  90 

93  87 

90  84 

87  81 

x=10 

96  87 

93  84 

90  81 

87  81 

84  78 

X=:ll 

90  81 

87  81 

84  78 

84  75 

81  75 

x^l2 

87  81 

84  78 

81  75 

81  75 

78  72 

X=:13 

84  75 

81  72 

78  69 

78  69 

75  69 

x=14 

78  69 

75  69 

72  66 

69  60 

69  60 

x=15 

72  66 

72  66 

69  63 

66  63 

66  60 

x=16 

69  63 

69  63 

66  60 

66  60 

63  60 

x=17 

66  60 

66  60 

63  60 

63  57 

63  57 

z^l8 

63  57 

63  57 

63  57 

60  57 

60  57 

x=19 

60  57 

60  57 

60  57 

60  57 

60  57 

x=20 

63  60 

63  60 

60  57 

60  57 

60  57 
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occccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


y=ll 

y=12 

y=13 

y=14 

7=15 

78  69 

75  69 

75  69 

75  69 

72  66 

z=2 

81  75 

81  75 

81  75 

78  72 

75  69 

x=s3 

87  81 

84  78 

81  78 

78  72 

75  69 

x=4 

90  84 

84  78 

84  78 

78  75 

75  69 

z=5 

90  84 

87  81 

84  78 

78  75 

75  69 

x=6 

90  84 

87  81 

84  78 

78  72 

75  66 

x=7 

90  84 

84  75 

84  75 

78  69 

72  66 

x=8 

84  78 

81  78 

78  72 

75  66 

69  63 

x=9 

84  78 

81  75 

78  69 

75  66 

69  63 

x=10 

81  75 

78  72 

75  69 

72  66 

69  66 

x=ll 

78  75 

78  72 

75  69 

75  66 

72  69 

x=12 

78  69 

75  69 

72  66 

69  66 

69  66 

xslS 

75  66 

72  66 

69  63 

69  63 

69  63 

x=14 

69  63 

69  63 

66  60 

66  60 

66  60 

z=15 

66  60 

66  60 

63  60 

63  60 

63  60 

X=:16 

63  57 

63  57 

63  57 

63  57 

63  57 

»  17 

60  57 

60  57 

60  57 

60  57 

60  57 

x=18 

60  57 

60  57 

60  54 

60  54 

60  54 

z=19 

60  54 

57  54 

57  54 

60  54 

6C  57 

x=20 

60  54 

60  54 

60  54 

60  57 

60  57 
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7=16 

7=17 

7=18 

7=19 

7=20 

X=1 

69  63 

66  63 

63  60 

63  57 

60  57 

x=2 

69  63 

66  63 

63  60 

60  57 

60  57 

x=:3 

69  66 

66  63 

63  60 

60  57 

57  54 

x=4 

69  66 

66  63 

63  60 

63  57 

60  54 

x=5 

69  66 

66  60 

63  60 

63  60 

60  57 

x=6 

69  63 

66  60 

63  60 

63  60 

60  57 

x=7 

69  63 

63  60 

63  60 

63  60 

60  57 

x=8 

66  60 

63  60 

63  60 

63  57 

60  57 

X=:9 

66  60 

69  60 

66  60 

66  60 

63  57 

X=10 

69  66 

69  66 

69  63 

66  60 

63  60 

Xsll 

72  69 

72  66 

69  63 

66  60 

63  60 

x=12 

72  66 

72  66 

69  63 

66  60 

63  60 

xsl3 

69  63 

69  63 

66  60 

63  60 

63  60 

x=14 

66  60 

66  60 

63  60 

63  57 

63  57 

z=15 

63  60 

63  60 

63  57 

63  57 

60  54 

x=16 

60  57 

63  57 

63  57 

60  57 

60  54 

x=17 

60  57 

60  57 

60  57 

60  57 

60  54 

x=18 

60  57 

60  57 

60  57 

60  57 

60  54 

x=19 

60  57 

60  57 

60  57 

60  57 

57  54 

x=:20 

60  54 

60  54 

60  54 

60  51 

57  51 

B.8 


The  data  which  follows  is  the  site  three  contoor  line  information  used  to  esti¬ 
mate  the  site  characterization  area  statistics. 


y^i 

y^2 

y=3 

y=4 

7=5 

X=1 

39  36 

39  36 

39  36 

39  36 

39  36 

x=2 

39  36 

39  36 

39  36 

39  36 

39  36 

x=3 

39  36 

39  36 

39  36 

39  36 

39  36 

x=4 

39  36 

39  36 

39  36 

39  36 

39  36 

x=5 

39  36 

39  36 

39  36 

39  36 

39  36 

x=:6 

39  36 

39  36 

39  36 

39  36 

39  36 

r=7 

42  36 

42  36 

42  36 

42  39 

42  39 

x=8 

45  39 

45  39 

45  39 

45  39 

45  39 

XssO 

45  39 

48  41 

51  42 

54  45 

54  45 

x=10 

48  42 

51  45 

57  48 

57  51 

60  51 

x=ll 

51  45 

54  45 

57  51 

60  54 

63  57 

x=12 

51  48 

54  48 

60  51 

66  57 

66  60 

x=:13 

54  48 

57  48 

60  54 

66  57 

69  63 

x=14 

54  48 

57  51 

60  54 

66  57 

69  60 

x=15 

51  45 

54  48 

60  51 

60  54 

63  57 

x=16 

51  45 

51  45 

54  48 

57  48 

60  51 

x=17 

48  42 

48  45 

51  45 

51  45 

54  45 

x=18 

48  42 

48  45 

48  45 

48  42 

48  42 

x=19 

45  39 

45  39 

45  42 

45  39 

45  39 

x=20 

42  39 

42  39 

42  39 

42  39 

42  39 
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y=6 

y=7 

y=8 

y=9 

7=10 

X=1 

42  36 

42  36 

39  36 

39  36 

39  36 

x=2 

39  36 

39  36 

39  36 

39  36 

39  36 

x=3 

39  36 

39  36 

39  36 

39  36 

39  36 

x=4 

39  36 

39  36 

39  36 

39  36 

39  36 

x=5 

39  36 

39  36 

39  36 

39  36 

39  36 

x=6 

39  36 

39  36 

39  36 

39  36 

39  36 

x=7 

42  39 

42  39 

42  36 

42  36 

42  39 

x=8 

45  39 

45  39 

45  39 

45  39 

45  39 

z=9 

54  45 

54  45 

51  42 

51  42 

48  42 

z=10 

60  51 

57  51 

57  48 

57  48 

54  42 

x=U 

63  57 

66  57 

63  57 

63  54 

60  48 

X=:12 

66  63 

66  63 

69  63 

69  60 

66  54 

x=13 

69  63 

69  63 

72  66 

72  66 

69  60 

x=14 

69  63 

69  63 

72  63 

72  66 

72  63 

x=15 

66  60 

66  60 

66  63 

69  63 

69  60 

x=16 

63  54 

66  54 

66  57 

66  57 

66  57 

x=17 

54  48 

60  48 

63  48 

60  51 

60  51 

x=18 

51  45 

51  45 

51  45 

5145 

54  45 

x=19 

48  39 

48  42 

48  42 

48  42 

48  42 

x=20 

42  39 

45  39 

45  39 

45  39 

45  39 

B-10 
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y=ll 

y=12 

y=13 

y=14 

y=15 

3P=1 

39  36 

39  36 

39  36 

39  36 

39  36 

x=2 

39  36 

39  36 

39  36 

39  36 

39  36 

X=:3 

39  36 

39  36 

39  36 

39  36 

39  36 

X=:4 

39  36 

39  36 

39  36 

39  36 

39  38 

x=5 

39  36 

39  36 

39  36 

39  36 

39  36 

x=6 

39  36 

39  36 

39  36 

39  36 

39  36 

x=7 

42  39 

42  39 

42  39 

42  39 

42  36 

z=8 

42  39 

42  39 

42  39 

42  39 

42  39 

x=9 

45  39 

45  39 

42  39 

42  39 

42  39 

z=10 

48  39 

45  39 

42  39 

42  39 

42  39 

x=ll 

57  45 

51  39 

45  39 

42  39 

42  39 

xsl2 

63  48 

57  45 

51  39 

45  39 

42  39 

x=13 

66  54 

60  48 

5142 

48  39 

45  39 

X=:14 

66  57 

63  54 

57  48 

51  42 

45  39 

x=15 

66  57 

63  54 

57  45 

48  42 

45  39 

x=16 

63  51 

60  48 

54  42 

54  42 

45  39 

x=17 

57  48 

54  45 

48  42 

45  39 

42  39 

x;=18 

54  45 

51  39 

45  42 

45  39 

42  39 

x=19 

48  39 

45  39 

45  39 

42  39 

42  39 

x==20 

42  39 

42  39 

42  39 

42  39 

42  39 
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y=16 

r=n 

7=18 

7=19 

7=20 

X=1 

39  36 

39  36 

39  36 

39  36 

39  36 

z=2 

39  36 

39  36 

39  36 

39  36 

39  36 

x=3 

39  36 

39  36 

39  36 

39  36 

39  36 

X=:4 

39  36 

39  36 

39  36 

39  36 

39  36 

Xs:5 

39  36 

39  36 

39  36 

45  36 

45  36 

x=6 

39  36 

39  36 

42  36 

45  39 

48  42 

x^7 

42  36 

39  36 

42  36 

45  39 

45  42 

x=8 

42  36 

39  36 

42  36 

45  39 

51  45 

z^9 

42  39 

39  36 

42  36 

45  39 

51  45 

xsio 

42  36 

39  36 

42  36 

48  36 

51  42 

x=ll 

42  36 

39  36 

39  36 

45  36 

48  39 

z^l2 

42  36 

39  36 

42  36 

42  36 

45  39 

x=13 

42  36 

39  36 

39  36 

39  36 

42  36 

x=rl4 

42  39 

42  39 

42  39 

39  36 

39  36 

xsIS 

42  39 

42  39 

42  39 

42  39 

39  36 

x=16 

42  39 

42  39 

42  39 

42  36 

39  36 

x=:17 

42  39 

42  39 

42  39 

42  36 

39  36 

x*18 

42  39  — 

42  39 

42  36 

39  36 

39  36 

x=19 

42  39 

42  39 

42  36 

39  36 

39  36 

x=:20 

42  39 

42  39 

42  36 

39  36 

39  36 
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